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2. 

1.  Introduction. 

This  paper  offers  a  nev;  approach  to  the  problera  of  dcternininc;  the  field 
propagated  in  a  non-homogeneous  atmosphere  by  an  arbitrary  source  of  eluctromacnetic 
radiation.  The  principal  problem  in  connection  v/ith  propagation  of  ultra-hi^h 
frequency  radiation  is  to  determine  the  refractive  effect  arising  from  inhomogen- 
eities  in  the  atmosphere.  Refraction  not  only  distorts  the  radiation  pattern  of  the 
source  but  produces  singular  or  anomalous  effects  such  as  propagation  beyond  the 
optical  horizon  and  trapping  of  radiation  at  definite  levels  so  that  the  intensity 
at  a  distant  point  vathin  such  d   level  is  far  gro.'ter  than  \/ould  exist  from  the 
normal  spreading  of  radiation  in  a  homogeneous  atmosphere.   Even  the  normal  atmos- 
phere produces  the  effect  of  "bending"  radiation  dov;n  beyond  the  optical  horizon. 
However , the  effect  of  a.  normal  atmosphere,  in  v/hich  the  index  of  refraction  de- 
creases lihearly  vdth  height  above  the  earth's  surface,  am   be  taken   into  account 
by  increasing  the  earth's  radius  by  a  factor  of  U/3  and  then  treating  the  atmosphere 
as  homogeneous  .   On  the  other  hand  the  dep/irtures  from  a  normal  atmosphere  pro- 
duced by  changes  in  temperature,  pressure,  and  './atcr-vapor  content, uhich  seriously 
affect  ultra-high  frequency  radio  wave  propagation , are  not  readily  evaluated.  The 

existence  and  magnitude  of  these  abnormal  effects  have  no;/  been  noted  so  often  ox- 

2 
perimentally  that  the  problem  of  deternining  theoretically  the  field  produced  by 

the  inhomogeneities  of  the  atmosphere  has  become  an  important  one. 

This  paper  apjjroaches  the  theoretical  problem  in  the  follov/ing  manner.   It 
aPGumes,  as  do  other  investigations,  that  the  cl^-ctromagnetic  parameters,  £  ,  yx  , 
and  (>    vary  only  uith  height  above  the  earth  but  not  in-  any  direction  parallel  to 
the  surface  of  the  earth.   It  acsuj-ies  farther  that  the  earth  is  flat.  (Correction 
for  this  ass-ujnption  and  the  accuracy  of  this  correction  factor  njre  discussed  in 
article  2.)   Under  these  assumptions  the  paper  shov/s,  first  that  solutions  of 
Maxwell's  equations  can  be  obtained  through  u/o  scalar  functions  V  and  'J  v/hich  satis- 
fy certain  partial  differential  equations  of  the  second  order  >[_2quat ions  (5. 2)  J  , 
f^^A^'he  ,f\inctions  V  and  \l   are  essentially  the  analogues  in  a  stratified  medi\im  of  the 
'Hvtfctz  vector  potential  './nich  serves  so  effectively  in  media  \/ith  constant  £  ,  //  , 
end  (T  , 


'Schelling,  J.C.,  Burrows,  C.R.,  and  Ferrel,  2.3.:.  Jltra-short  ..'ave  Propagation, 
Proc,  I.H.Z, ,  V.  21,  1S33,  PP  1^27-^63.   Also  2ckersley,  T.L. :  Ultra-short  -/ave  He- 
fraction  and  Diffraction,  Jour.  Jnst.  ilec.  2ng.  ,  V.  '60,    1937,  PP.  236-30U. 

"'i'iany  of  the  references  cited  later  in  other  connections  give  experimental  data. 
Systematic  experimentation  is  being  done  also  by  the  Georgia  School  of  Technology. 


vi'hen  the  source  is   a  plane  uave,  expression;  for  V  and  V/  v;hich  satisfy  the 
partial  differential  equations  are,   at   least   theoretically,   readily  obtained.     To 
treat  next  any  arbitrary  source  of  electromagnetic  radiation  it   is   supposed  that 
the  source  can  be  expressed  in  a  homogeneous  mediun  as  a  superposition  of  ordinary 
plane  uaves^    |_Sxpression   (5.7)  >     •     This   supposition  is   certainly  correct  for  the 
most  important   single  type  of   source,   namely  the  electric  dipole,    and  for  many 
others.     By  a  mathematical  argument  v;hich  is   fundamentally  simple  it   is   then  sho;m 
that   the  radiation  field  produced  in  the  non-hcmogeneous  medium  by  an  arbitrary 
source  can  be  obtained  by  a  superposition  of  the  radiation  fields  produced  by  the 
plane  v;aves  v/hich  comprise  the  source.      The  amplitudes  v/ith  which  the   individual 
plane  uaves   enter  into  the  mathematical  representation  of  the  field  of  the  arbitrary 
source  when  in  a  homogeneous  mediina.of  course  enter  into   the  determination  of  the 
field  in  the  non-homogeneous  medium.     That   is   to  say,    the  property  v;hich  distinguishes 
one   source  from  another  resides   in  an  amplitude  function  v.'hich  depends  upon  the 
directions   of  the  constituent  plane  v/aves.. 

This   argument   is  applied  first   to   the  case  \/here  a  source  is   located  be- 
lov;  a  semi-infinite  non-homogeneous  mediiim  and  the  field  in  the  medium   is  then 
determined.     The  more  general  case  of  an  arbitrary  source  located  bet\;een  t\;o  semi- 
infinite  media  is   treated  next.     This  section   can  be  specialized  to  treat   the  case 
of  a  source   of  radiation  located  above  the   earth.     The  raain  result   of  the  paper  is 
embodied  in  formulas   (7.S)   and  (7.9),   vhich  give  the  expressions  for  V   and  ¥, 
These  expressions,   taken  in  conjunction  v;ith  lormalas    (5.3)»  permit   determination 
of  the  six  components    of  the    electrouagnetic   field  from  V  and  \I,     A  number  of  special 
cases  are  './orlced  out   to  illustrate  the   theory^    although  they  do  not   necessarily 
represent  any  actual  atmospheric  conditions. 

It  v;ill  be   seen  that  the  main  result  ,1  Equations      (7.S)  and  (7.9)  J, is   in 
the  form  of  a  double  iritegral  of  a  complex  function.      The  problem  of  evaluating 
this  integral  for  any  particular  atmospheric  lav  of  variation  of  index   of  refraction 
vfith  height  can  be  a  difficult   one.      Some   of   the   examples   of  the  paper  illustrate 
this  point.     However,    difficulties   in  the    evaluation  of  the   integral   are  to  be  ex- 
pected.    As   our  discussion  below  will   indica.te,   the   one   satisfactoiy  approach  which 
has   thus  far  been  made  to  the  problem  of   anomalous  propagation  has  produced  tremendous 
calculation  problems.      The  problem  of  evaluation  of  the    integral  laay  be  considerably 
simplified  by  the  use   of   reasonable  assumptions.     For   example,    it   is  practical   to 
use  the  fact   that   the  variation  of  the   index   of  refraction  \n.th  wave   length  is  very 
small.     As  a  matter  of  fact   this   assumption  is  i-iade  by  other  investigators  at  the 
very  outset  of   their  v/ork. 


The  present  approach  to  the  problem  of  anomalous  propagation  is  an  alterna- 
tive to   the  aoproach  initiated  and  developed  by  several  uorkers   in  the   field,   not- 
a"bly  Furry,   Hartree,  Pearcey,   Pryce  ?,-id  Pclceris''   .      The  latter  approach  may  be 
described  briefly  as  follous.      Under  the   assiir.iption  that   the  variation  in  the  index 
of   refraction,    n,    is   snail  compared  to   the  './ave  length,    A',    the   steady  state   dlectro- 
iiiagnetic  field  can  be   obtained  froi-i  the  function  \i/(x,y,z)  v;hich  satisfies  the  \/ave 
eqxiation 

S?^VJ  +  k^  n^Y     =  0  (1.1) 

uhere  k  =  -^  ,  Jor  a  vertical  dipole  y   represents  the  vertical  conponent  of  the 

electrical  Hertzian  potential  f-onction  of  the  dipole,  and  for  a  horizontal  dipole 
y  represents  the  vertical  conponent  of  the  .Tiagnetic  Hertzian  potential  function  of 
the  dipole.   (It  can  be  sho\;n  that  the  field  intensity  |El  is  proportional  to  [V/j 
in  each  case, ) 

This  equation  must  be  solved  for  ^t/  subject  to  several  boundary  conditions 
v/nich  arise  froa  the  fact  that  the  source  is  a  dipole,  that  the  ground  is  a  perfect 
conductor,  and  that  the  \;ave  i.Tust  be  outgoing  at  infinity.   It  is  then  assumed  that 
the  earth  is  flat,  correction  being  r^de  for  this  assumption  by  the  use  of  a  modi- 
fied index  of  refraction  (see  Art.  2). 

The  solution  of  equation  (1.1)  is  obtained  in  the  fona  of  an  infinite 
series 

y   =     -  i^  \_  H^'  i\-)  \i\)   Ujh^)  (U2) 

where  d  is  horizontal   distance  from  the  source  and  h     and  h     are  the  heights  of  the 

'■>':t,Tansiaitter  and  receiver  above  ground.      The  functions  U     and  the  quantities  k 
•'■Ma,  I  u  T  m  ^  m 

must  satisfy  the  differential   eouation 

Mk2K2_k2^  U  =0.  (1.3) 


dh^  ^    ° 


-./here  IT  is  nov;  the  modified  index  of  refraction  as  a  function  of  height  h  above  a 
flat  earth  and  U^  must  satisfy  some  boundary  conditions  arising  out  of  the  condi- 
tions on  U'  .  Because  of  these  boundary  conditions  the  values  of  1:  are  limited  to 


m 


certain  characteristic  values.   Hence  only  certain  characteristic  functions  U  .each 

m' 
'..'ith  its  corresponding  k  ,  -nay  be  Used. 


The  essential  theory  of  this  approach  is  given  by  Furry  U.H..:  Theory  of  Character- 
istic Functions  in  Problems  of  Anomalous  Propagation,  Radiation  Laboratory  Report 
Tf620,  5^eb.  1945. 


5. 

Tais  alternati-co  method  calls  for  finding  the  characterictic  functions 
and  characteristic  values  of  the  second  order,  ordinary  differential  equation  (I.3) 
v;hich  cannot  "be  solved  exactly,  except  in  very  special  cases  wnich  do  not  happen  to 
represent  real  atmospheric  conditions.  The  determination  of  the  characteristic 
values  and  chara.cteristic  functions  is  a  major  problem  of  calculation  and  calls  for 
the  use  of  advanced  mathematical  methods.   Because  of  the  complexity  of  the  problem 
many  separate  methods  have  been  investigated  and  compared.  Among  the  methods  tried 
there  are  the  phase  integral,  or  3.',/.K.  method;   the  perturbation  method  ,v;hich  has 
been  found  to  be  applicable  to  complex  eigenvalues  in  v/hich  the  imaginary  part  is 
large  (this  corresponds  physically  to  the  case  vhere  the  energy  represented  by  the 
terns  of  the  above  series  leaks  copiously  from  the  ducts  formed  in  the  atmosphere); 
Rayleigh's  method,  v/hich  appears  x,o   be  applicable  best  to  real  eigenval\;igs  (physical- 
ly this  corres]ponds  to  trapped  modes);   the  variational  method /mich  seems  to  work 
v;ell  for  real  esid   couiplex  eigenvalues  uhich  correspond  to  viaat  may  be  physically 
described  as  some\;hat  leaky  modes;   ctnd  finally,  the  differential  analyzer,  \/hich 

merely  resorts  to  a  machine  to  handle  tne  labor  and  reduce  the  time  reouired.   Some 

k 
discussion  ol  these  methods  v/ill  be  found  in  paijers  by  G.G.  t'iacfarlane  and  C.  L. 

5 

Pekeris  • 

It  is  fair  to  say  that  all  these  methods  have  been  only  partially  success- 
ful in  that  any  one  method  works  for  crdy  very  special  la\/s  of  variation  of  x*  uith 
height  and  then  only  at  the  cost  of  tremendous  calculations..  For  some  la\/s  of  varia- 
tion no  methods  have  been  successful  in  getting  good  approximations  to  the  character- 
istic values.   Indeed, to  secure  approximations  1,0  i/hich  one  can  attach  some  degree 
of  assurance, it  appears  to  be  necessary,  since  the  correct  values  are  not  knoun,  to 
use  several  methods  of  calculation  simultaneously  aiid  to  rely  upon  agreement  of  the 
results,  if  srach  agreeiient  is  obtained. 

When  the  characteristic  values  and  characteristic  functions  are  obtained 
to  some  degree  of  approximation  they  must  be  substituted  in  equation  (1.2)  and  the 
series  there  summed.  Ihis  summation  is  often  feasible  to  a   useful  degree  of 
approximation  only  in  some  ranges  of  d  values. 

jj- — 

'kacf arlsne ,    G.G.  :  A  Variational  liethod  for  Determining  Eigenvalues   of  the  V/ave 
Equation  Applied  to  Refraction.      Cambridge  Fhilosonhical   Society  Proc.,Y.   U3  Part   II, 
April,    19^7. PP. 213-219. 

•Pekeris,    C.L.;   V/ave  Theoretical  Interpretation  of  Propagation  of  10  Cm.   and  3  Cm. 
V/aves   in  Low-Level  Ocean  jDucts,   Proc,    I.H.3.,   V.    35,   May   I9U7,  pp.   U53-I162. 


6. 

The  theory  oi  this  paper  \ri.ll,  as  renarhed  above,  lik&i/ise  entail  calcu- 
lation problems.   At  the  present  tine,  all  that  can  be  said  is  that  it  has  been 
possible  in  other  proble;ns  of  applied  .aather-i-tics  to  obtain  good  approximations  to 
the  values  of  integrals  by  famili?.r  i/iCthcds,  such  as  the  niethod  of  steepest  descent  , 
\;here  series  solutions  have  led  to  tremendously  lon^  calculations  and  to  very  slov; 
convergence. 

The  theory  of  this  paper  is  considerably  more  general  than  other  approaches 
in  that  the  radiation  source  is  arbitrary  cC.d.   the  medium  belov;  tlie  source  is  arbit- 
rary.  Other  investigations  have  been  limited  to  the  dipole  source  for  the  reason 
that  the  propagated  field  is  generally  defired  at  points  near  the  ground  and  at 
large  distances  from  the  transmitter..  It  is  Ir.o'.m   Tihat  radiation  in  non-homogeneous 
.atmospheres  reaching  distant  points  near  the  gro.md  emanates  from  the  source  in 
directions  deviating  less  than  a  degree  from  tn^  horizontal.   Hence  it  is  argued 
that  by  computing  for  tne  dipole  the  ratio  of  energy  per  unit  area  received  to 
energy  per  unit  area  passing  horizontally  at  some  fixed  distance  from  the  trans- 
mitter, one  could  compute  the  field  received  from  any  trarismitter  by  measuring  or 
calculating  the  energy  produced  at  the  fiied  distrr.ce.   However,  there  are  many 
applications  of  a  theory  of  propagation  to  v;hich  this  argument  does  not  apply.   For 
example,  to  obtain  the  coverage  pattern  of  a  radar  station  for  purposes  of  airplane 
detection,  communications,  and  landing  systems,  requires  more  than  a  knc\/ledge  of 
the  field  propagated  along  the  surface  of  the  earthy  That  the  coverage  pattern  of 
a  radar  set  can  be  considera,bly  modified  by  anomalous  atmospheric  conditions,  has 
been  stressed  by  Freehaf er. 

The  remarks  of  the  last  f.jv;  paragraphic  indicate  reasoriS  for  undertalcing  a 
nev;  a'O'oroach  rather  than  for  •oursuing  further  an  apriroach  \;hich  has  been  somewhat  ex- 
\'ph:T&d   already. 

It  should  be  noted  tha.t  a  theory  of  propagation  through  stratified  media 
has  application  to  the  design  of  interference  films  (as  on  glass)  so  as  to  eliminate 
reflection  or  to  eliminate  elliptical  polarization  from  totally  reflecting  surfaces. 

'See,  for  example,  Stratton,  J.  A,:  electromagnetic  Theory,  i;cGra\/  Hill  Book  Co., 


Inc.,  xl.  Y.,  lQl+1,  p.  368. 

?reehafer,  John  3. :   The 
Hadiation  Laboratory  Report  =if'^'''tf. 


?reehafer,  John  3.:   The  ^iffect  of  Atmospheric  Refraction  on  Short  Radio  \Iavez, 


7. 

2.  The  Use  of  a  Flat  Sarth 

In  this  paper,  as  in  other  investigp.tious  of  the  problem  of  propagption 
of  high  frequency  v/aves  in  non-homogeneous  a tnio spheres,  it  is  supposed  that  the  sur- 
face of  the  earth  is  flat.   The  errors  introduced  thorehy  are  corrected  somewhat  "by 
a  device  now  standard  Icio'i/n  as  the  modified  index  of  refraction.  As  remarked  in 
article  l,it  is  assuned  conmonly  that  the  index  of  refraction  varies  only  with 
height  above  the  surface  of  the  earth.   One  then  replaces  the  variation  of  the  in- 
dex of  refraction  with  height  above  the  surface  of  tlie  spherical  earth  by  a  aodi- 
fied  index  of  refraction  which  is  given  by  the  formula 


'•<')  =  ^^  (=-^) 


where  n(r)   is   the  actioal  variation  of  the    index  n  v/ith  r,    the  distance   from  the 
earth's  center,  -and  b   is   any  fixed  value  of  r  usioally  taken  to  be     a,    the   earth's 
radius.      If  v;e   replace  r  by  a  -♦•  z  in  for^rala   (2.1),    since  a  is   constant,   we  have  a 
variation  of   index  of  refraction  v;ith  z  alone,    z  now  being  the  height    above  the 
earth's   surface   in  a  rectangular  coordinate   system.      In  this   coordinate   system  x 
and  y  represent    distances  along  the  flat   surface.     Any  point   on  this    surface  corres- 
ponds to   the  point  on  earth  having  the   same  direction  and  the  same  distance   (along 
the  great   circle  route  on  the   earth's  surface)   from  the   origin. 

Physically  expressed,    tiie   justification  for  the  modified  index  of  re- 
fraction  is   roughly      that  .instead  of  having   the   earth  curve  av;ay  from   the   radiation, 
the  index  is   changed  in  such  a  way  as   to  nake  the  rays   curve  av;ay  from  the  earth. 
The  mathematical  justification  for  the  use   of  this  modified  index  of  refraction, 
v/hich  corrects    only  approximptely  for  the    "flattening"  of  the   earth,    is  given  most 
clearly  in  a  paper  by  Pekeris    ,    wherein   estimates   of  the   validity  of  the   approxima- 
tion are  also  given.      Essentially  the  approximation  is  valid  out    to   ranges  of  the 
order  of   the   radius   ox    the   earth  , for  wave   lengths   larger    than  a  centimeter, and  to 
heights   of  the   order   of  1,000   feet.      Ii   the   rarige   or  height   is   exceeded  or  the  v;ave 
length  is   decreased  the   approximation  becomes  p-jver.      Actual   data  are   given  by 

Pekeris    in  his    report.      The   earth-flattening  a.pDroximation   is   considered  for  -our- 

9 
■ooses   of  ray-tracing  by  j'reeimfer.'^ 


'Pekeris,  C.    L. :   Accuracy  of  the  jiarth-Plattening  Approximation,    Col.    U.   Math,  Phys. 
Group  Report  ^3,   Apr.    15,    I9U6.     Also  in  Phys.    Rev.   V.    70,    ISkh,   pp.    518-522, 

9  _ 

/'jreehafer,  John  '£.:    loc.cit.,  see  page  5  ^^'^  particular. 


8. 

It  is  also  necessary  in  the  case  where  the  earth  is  treated  ?s  an  imper- 
fect conductor  to  correct  for  the  index  of  refraction  in  the  earth.   ?or  a  conductor 
the  index  of  refraction  n  is  given  hy  ^T^  IT  \;here  c     ~         -     — ,  •  .and  this  n  nust 

he  modified  "by  altering  £   and  S   to  produce  relation  (2.1).  Ko  correction  is  nade 
foryu-jfor  it  is  practically  one  in  all  atmospheres  near  the  surface  of  the  earth  and 
for  short  distances  in  the  earth. 

The  theory  of  this  paper  and  the  illustrative  examples  deal  directly  with 
i   ,  jU  ,    and<y  rather  than  v;ith  the  index  of  refraction.   iTo  mention  is  made  as  to 
v/hether  these  values  are  modified  or  not.   In  the  application  of  the  theory  to 
anomalous  propagation  over  a  flat  earth,  £.  and  <^  would  "be  modified,  i/hereas  in  the 
application  to  a  prohlen  of  interference  la.yers  they  \;ould  not  he.   In  the  former 
case  some  of  the  illust^-ationc  of  the  theory  given  in  this  pa.per  must  he  properly 

interpreted.   For  example  if  the  medium  used  for  imrposes  of  illustration  is  sup- 
posed to  consist  of  plane  parallel  layers  v/ithin  v/hich  c  ,  yU  ,  and  C  are  constant, 
this  situation  over  a  flat  earth  corresponds  to  one  in  which  the  actual  variation 
of  the  index  of  refraction  with  height  ahove  a  spherical  earth  is  of  the  form  — 

v/ith  k  different  from  layer  to  la^'cr.  Just  which  cases  of  actual  variation  of  n 
with  height  above  a  spherical  earth  :aay  he  usefully  approximated  by  a  series  of 
functions  of  the  form  —  is,  for  the  moment,  irrelevgmt.   In  fact,  the  theory  does 
not  require  that  an  actual  atmosphere  he  treated  by  approximating  it  with  a  finite 
nvui'iber  of  layers  of  definite  thickness. 

Beicause  the  approximation  introduced  by  the  modified  index  of  refraction 
is  not  good  for  \/ave  lengths  less  than  1  cm  or  :^or  heights  above  about  1,000  feet,, 
it  seems  advisable  for  these  ranges  to  develop  a  theory  for  a  spherical  earth  directly 
and  avoid  this  earth  flattening  approximation  entirely.   Some  throught  has  already 
been  given  to  a  theory  of  propagation  for  a  stratified  mediun  above  a  spherical 


earth.  Stratification  in  this  case  ;.ieans  that  L  ,  jJL  ,   and  CT  are  functions  of  r 

/~2    2   ""2 

only  where  r  =  \/x  +  y  +  z 
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7     Some  Theory  Concerning  Majcv/ell's  Sqiia.tions. 
The   electromagnetic  vectoro 

H=    (H^.    H^.   H3) 
satisfy  waxv/ell' s   eqiiations 


-£-3^   -  curl  r-  =     E  /^  ,  % 

c      t  c  (j;^l) 


V^H,   +  curl  S  =  0 
c      t 


Under  the  influence  of  a  periodic  Eo-orce  a  field  S,  H  is  established 
sriodic  in  time  and  can  he  represer 

'"^^*'*'         •  -    £  =  Heal  Part  of  u  c~^"*^ 


Avhich  is  periodic  in  time  and  can  he  represented  in  the  form 


-ict  (3.2) 

H  =  Real  Part  of  v  e 


\/here  the  complex  vectors  u  and  v  satisfy  the  equations 
„■•  curl  V  +  i  k  £  u  =  0 


(3.5) 

curl  u  -  i  k  ^v  V  =  0 


where  k  =  -^=  -r^  and  t  =  c  -  1 


Me   shall  assume  in  the  follov;ing  that  6=0   and  thus  discuss  only  the 


equations 


^  S  -  curl  H  =  0 


■A  E,  +  curl  2  =  0 
c   t 


(3.U) 


'ie  note,  ho\;ever,  that  from  a  tiue-periodic  solution  S,  H  of  these  equa- 
tions a  correct  periodic  solution  of  equation  (3.I)  '..'ith  6"   t  0   can  be  found  simply 

*        U~  '^ 
by  replacing  the  real  quantity  £  by  the  con-olex  quantity  i.     =  £.-   '     i. 

At  any  point  x,y,z  \/here  there  is  no  source  the  electromagnetic  vectors 
E  and  H  satisfy  the  additional  equations,  div  £  H  =  div  nE  =  0,    so  that  the  complete 
sy&tem  of  equations  is 

a)  -^E^  -  curl  H  =  0,  c)  div  £  E  =  0 

V  U  0.5) 

"b)  -^H^  +  curl  2  =  0,  d)  div  yi<  H  =  0 
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3y  using  thr;   fnct    that,   the   divergence  ol"   tlie   curl   oi"  any  vector   function 
ic   0,    it    is   seen   that   the   second  set   of   equations    is   a  consequence   of  the   first   if 
div  c  li  -   div   u 'd  =  0  at    ar^  pm'ticulnr   ti..iQ    z.      In  case   of   time-periodic   solutions, 
the   latter   t\;o    equations  are   a  direct   consequence  of    (3o)- 

..''e   can  replace   the    equations    (3.5)  "'^y  second-order  equations    for   the  com- 
ponents  of   the  vector  S  or  of   the  vactor  H.      3y  talking  the   curl   of  equation  'o)  and 
substituting   it    in  the   tine   derivative   of   equation  a)   \/e   obtain 

-^'4^    Z^,  +     Acurl     -     curl  2  =  0 
2        tt  /x 

c  ' 

Similarly ,\;o  obtain  .  (3.. 6) 

— ^    r-  ,   ■*-    c  curl      ;.-     curl  h  =   0 

(d  X/  '-J  i^ 

c 

If  a  solution  Z  of   the   first    eouation  is  hnov/n  for  I'hich  div    £.  3  =  0  ,then 

C 

the  associated  magnetic  vector  H  ca.n  be  found  by  /^K     -  -^curl  S.      Similarly , if  a 

solution  K  of  the   second  eoixation   is   :-:nov/n  for  uhicli  div    /lR  =  Othen  3   follous   from 


r 


c 

c3.    -^curl   H. 

IP 
V/ith  the  aid  of   the   vector   identities 


ucurl  —     curl  ii  =  -   ZiZ  ■*■  grad  div  3  +   /.^(grad  —  )  x  curl  E 
and  ^  div    6,     3  =  0  =   div  3  --  —  (grad  £    )    .    3 

ue  can  \/rite  equations    (3.S)    in  the  form 

^^  2++     -  £i3     =  grad  (p   •    3)     +  q  x  curl  3 
c  \ 

^Ay  H^   -  A  H  =  grad  (q  •  H)  +  p  x  curl  H 
c 


(3.7) 


(3.3) 


where  the  vectors  p  and  o^  are  defined  by 

p  =  grad  (log  c  ) 

q  =  grad  (log  ;^.  ) 

Hote   that   the   right   sides   of   eqiJiation    (3.7)   contain  only  derivatives   of   first   order. 
In  general,    the   six   eq-'ii.tions    (3.7)    involve   all  coniponcnts    in  each  equa- 
tion.     The   simplication  in   the  case  of  a  stratified  medium   is   that   t\;o   equations   can 
be  obtained,    one   involving  only  one  component   of  3  ,and  the   other     only  one  component 
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See  Stratton,  loc.  cit. ,  p.  60U,  relations  (9)  and  (10). 
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of  K.  These  components  are  the  z-coinponents  E-,  and  r'  in  the  case  of  a  stratified 
nediun.  This  is  shoun  in  the  next  section. 

1|.  Representation  of  3  and  H  in  a  Stratified  Medina. 

\!e   now  specialize  the  results  of  article  3  'to  the  case  of  a  stratified 
medium.   3y  such  a  medium  ve   mean  one  in  uhich  c  ,  /j.  ,  and  O" ,   vary  in  one  dimension 
only.   Since  this  paper  replaces  the  spherical  earth  "oy  a  flat  one  ty  making  the 
customary  change  in  index  of  refraction,  the  direction  in  ;;hich  £,//.,  and  cf  vary 
will  be  taken  to  "be  perpendicular  to  the  flat  earth  and  this  direction  indicated  "by 
the  z-coordina.te  of  our  rectangular  coordinate  system.   It  mast  no\/  be  understood, 
however, that  ^  , /^  ,  and  QTare  the  modified  characteristics  of  the  atmosphere  in 
accordance  uith  the  discussion  of  article  2. 

Vi'e  have  in  this  case, from  equation  (3.0), 

P=4-k;      ^=#- 

where  primes   denote   differentiation  \/ith   respect   to    z,    and  k  is  the  unit  vector   in 

the   z-direction. 

Hence,  from  eq\iations    (3.7)f 


J±    3        _  ^3  =  grad   (4-  EJ  +  ^  (k  X  curl  Z) 


H^^   -  £>  H  =  gra.d   (-'-^  H, )   -f-  -■?-  (k  x  curl  H) 


(U.l) 


C 


In  particular,   by  equating   the  third  components   of  these   tv/o  vector   equations,   v;e 
find  that  S,  and  H^  satisfy,    respectively, 

i^  ■?      _  A  7    =  JL  (J-     V  'l 

c^   3^t  3       az    ^    -5  (^^2) 

^   3tt     ^''3         az  ^  A  3 

Our  next  problem  is  to  determine  the  remaining  components  E, ,E„,  K  and  H, 

-I.      2       1  c 

?or  this  purpose  v;e  use   the   equations  div    i  E  =  0,    div  yUH  =  0  and  the   relations   ob- 
tained for  the  components  E^  and  E^  from  the   first  two   of  equations    (3.5), 

These  four  relations  are: 


5-<-         ciy  £  Sz  '     Sx  ay      ~  ju        dz 

_^  _  _2    =  ^    __!_  .     ,,  1  _    2  ^  _  £    i_ 

ay       d'j^  c       o"t  '      c>y  c^  c      dt 


(^.3) 
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Solutians  of  these  eo^uations    can  "be  foiin.d  as   follo\;s,     V/e   assume  that  3,  and  H^  can 
le  reDresented  "by  tv;o   functions  V(x,y,z,t)  a.-:d  './(x,y  ,z,t )  ,as  follov;s: 

£E,  =  -  (If;    .    Ifl) 

(U.U) 


^  dx  By 


This  assumption  does  not  introduce  any  loss  of  generality,  since  to  any  given  pair 
2  ,  H,  many  functions  V,  './  can  "be  constructed  so  that  equation  C^.  U)  is  satisfied. 

3y  introducinii;  equations  (U.U)  in  equations  (U.2)  v/e  find  that  if  V  and  './ 
satisfy  the  equations 

.   1^  V,,  =  V    +  Y    +  6(4  V  ) 
2   tt    XX     yy      t  z 

""  (U.5) 

^2  tt    :ac    yy    ^^'^^  z  2 

then  the  values   of  S     and  H     given  hy  relations    (U.U)  uill    satisfy  equations    (U,2). 

^j'ith  the  representation  (U, U)   of  2^  and  H,  we  now  get    solutions   of  e'-uations 

<SK    =  V       -  ^  V/   ^  oE-,    =  V;'       +/^V  ^ 

T.         X2       c     yt  /^   1         xz     /c     yt 

(1+.6) 

^    2        yz       c     xt  /^    2        yz     'c     xt 

as  one  easily  verifies  "by  direct  substitution. 

On  the   assuTiiption   tliat   i-iaxv/ell's   equations   hold  v/e  v/ere   led  to   equations 
(^.^),    (^.5)»   ^^^   (U.6).      For  the  piirpose  of  this  paper  v/e  must  show  next  that    the 
two  vectors  E  and  H  defined  hy  equations    (U.U)   and   (U.6)   satisfy  Maxv;ell's   eq-uations; 
that   is,v;e  must    shov;  that  v/e   can  reverse   our   steps. 

Theorem.      If  V   and  ',J  are   solutions   of   the    follouing  partial   differential   equations, 
Vv'herein  £  and  m  are  functions    of  z  only. 

^       tt         XX  yy        '-'c       z\ 

(U.7) 


(U.3)  ty 


2      tt  XX  yy        r    Li.     z 
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then  a  solution  E,   H  of  Kaxv/ell's   equations    is   ottc-ined  "by   the    foraulae: 

-     ----  ■■  —  ■■i  -       ■■  ...■  .  —    ■---      —    __ 


3  XX        yy       '  '^   j>  XX        yy 

Proof.  \ie   note  that  (U.S)  oay  1)6  written  in  vector  fonn 

£.S  =  curl  A 

LLo.  =  curl  5 
\/here  A  and  3  are  the  vectors 


(U.9) 


A  =  (V  ,  -  V  ,  -  ^  \I    ) 
y'    x'   c  t 


B  =  i\l^.    -  \>ft^t^ 


C4.10) 


Since  the  divergence  of  the  curl  of  any  vector  function  is  zero,,  it  is 
evident  that  div  £  S  =  div  i£  H  =  0. 

Purthermore,  from  equations  (U.2)  ana  (U^IO) 

^  S,  -  curl  H=  curl  (i  A^  -  H)  ^^^^^^ 

and  "by   equations    (U.8)   and   (U.IO) 

-  A^   -  H  =   (-  i  '.;      ,_  1  U    ^    ,-  -L  \;   ^   +  i  V/       *  1  W      ) 
c     t  ^     xz'  yu     y?         ^2     tt     /:<.     XX     ^    yy' 

=   (-  ,-  W      ,   --\I      ,    -   (-^  )„] 
V  /x    xz'     /^    yz'       '^i^   'V 

■fay  the   second  of  relations   (U.y); 

=  grad  ( -)    . 

Again  "by  relations   (U.8)  and  (U.IO) 

I3     +E=(—       ^-     /^V       -iv       -Iv      ) 
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5y  the  first  of  relations  (U. 7)  this 

=  grad  {-^)      . 
c 

Since  the  curl   of  the   gradient   of  any  vector  function  is    zero,    ue  have  by  equation 
(U.ll) 

^\ 

curl  H  =  0 

c 

/-  K^   +  curl  S  =  0 
c     t 

v/hich  proves  our  the  or  en. 

Our  above  consideration  does  not   shov/,    of  course,    that  every  solution, 
E,   H,    of  I'Iax\/ell's  equations  , for  which   div  £  3  =  div  it  K  =  0  and  \;herein  6-  and /-t 
are  functions  of  z  only,  can  be  represented  in  the  form  of  equation  (U.8).      Kov/ever, 
it    is   not   difficult    to    see  triat   this   is   the   case.      It   is  aade  plausible  by   the 
following  reasoning. 

The   solution  of  ;-;ac.\;ell's   equations    is  uniquely  determined  by  assign- 
ing initial  values   to    the  six  conoonenos  of  'jL  and  H  at   t  =  0.      However,    these 
initial  values  nust   satisfy  th-j   condition  div  £3  =     div  aj  H  =  0, which  reduces   the 
"degree     of  freedom"   in  essence  to   four  arbitrary  functions,      Exactly  four  functions, 
ho\/ever,    determine  the   solutions  V  and  \i  of  equ.-tion  (^.7)  namely  the  values   of 
V,   '.[,   V.  ,    and  \I.    at  t   =  0. ''■■'■ 
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Our  formulas  contain  as  a  special  case  f or  £  =  ix  =   1  and  VJ  f  0  the  Hertzian 

radiation  from  a  di-oole,    ■  r        ,  .j.  \  ^\ 

^i(r  -wt)  ^^.  9j 

Letting  V  =  we  find  from  eouations  (U.IO),  -fU.lH  and 

(4.1)  ^ 

(•^ri-b  that 

E  =  curl  (iV  -  jV  ) 
-  y   -  X 

H  =  —  curl  V,  k  =  curl  tt^  .  tt  =  -  Vk. 
c       t—        t      -c    — 

TT  is  the  Hertz  vector  of  a  dipole. 

That  the  field  of  a  more  general  Hertzian  vector  tt  can  be  represented 
by  eq-oation  t^+»-i2J.  follo\/s  from,  the  fact  tiiat  curl  grad^^'  =  0  and  that  to  a  given 
vector  TT  a  vector  tt'  =  tt  •♦■  grad<i?  can  be  constructed  such  that  tt'  has  the  special 
form  (Uy  .  -\   ,/f\). 


The  theoren  of  this  article  deserves  special  attention.   It  s.ays  thn.t 
for  stratified  uiedia,  i.e.,  nedia  in  './hich  ^,  M,  and    CJ  functions  of  z  only,  the 
solution  of  i'iax\;ell's  equations,  \;hich  involve  6  functions  with  at  least  tv;o  func- 
tions in  each  equation,  can  "be  reduced  to  the  proole^i  of  finding  two  functions, 
naxaely  S,  and  H, ,  with  a  separate  equation  for  each.   This  result  is  reached  in 

equations  (U.2).   Additional  theory  then  redxices  the  protlea  of  finding  S,  and  H, 

7  3     3 

to  that  of  solving  equations  (U.^),  one  for  V  and  one  for  './.   Formulas  (U.S)  then 

give  the  reniaining  four  functions. 

Standard  electrouagnetic  theory  yernits  the  reduction  of  the  problem 

of  solving  i4ax\;ell's  eq-'oations  to  the  soliition  of  a  systea  of  differential  equations 

each  involving  only  one  function  only  in  the  case  that  c. ,.  -ii  ,  and  (Tare  constant. 

Hence  the  result  of  this  x;aper  extends  standard  theory  at  least  to  stratified  media. 

5»  Solutions  for  V  and  V/  in  the  case  of  a  stratified  medium^ 

Ve  have  seen  in  tne  preceding  sections  tliat  the  electromagnetic  field 
in  a  stratified  medium  can  be  represented  "ith  the  aid  of  two  scalar  functions  V 
and  W.   These  functions  satisfy  the  partial  differential  equations 

%^t  =  ^x:c*^yy  -^  4^z), 

(U.7) 

^"2  'tt    -xx       -yy  /^^n\\ 

C  ' 

The  electromagnetic  vectors  E  and  H  then  are  oUtained  by  the  formulae  (U.S).   V/e 
proceed  to  obtain  functions  V  and  V/  \/hich  satisfy  these  equations,. 

We  do  not  assiane  in  the  follov;ing  that  c(z)  and  ju,{i)   are  continuous 
functions,  though  \ie  shall  require  that  they  are  sectionally  smooth.  At  any  surface 
of  discontinuity  the  functions  2-^,   3p,  ^S,  and  K,.,.  P",  /^  H_  must  be  continuous. 
ue   recognize  by  relations  (U. 8)  that  this  is  the  case  if 

r 

the  functions  V  and  'il   and  the  functions  :r  V  and  -17  are  continuous   (5.1) 
C  z    ^  z 

V/e  shall  therefore  seek  a  V  and  ;/  \;hich  satisfy  these  conditions.. 

Ue  shall  confine  ourselves  in  this  section  to  solutions  S  and  E  v/hich 
are  periodic  in  time  and  thus  have  the  form 

Rize-i^^l.  and   R  1  H  e-^^H 
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[■  where  E(x,y,z)  and  H(x,y,z)  are  complex  vectors   independent   of  t.     '.!e  obtain  such 

t  solutions  \n.th  the  aid  of  functions  V,   \1  uhic.i  have  the  form  V(x,y,z)e~'^  , 

.  VKx,y,z)e  ,   uhere  V  and  './  nov;  are  independent   of  t  and  satisfy  the  equations 

'  obtained  from   equation   (U.7)   "by   substituting;   these   ne\-r  forms    for   the   V  and  \!  there. 


(i  V   )       ■*-  k^n^     = 


V      +  V        +    £C^  V  )       +  k  n^    =  0 
XX         yy  ^     z   ,  ^^^^^ 

W       +  V/         ■»•  A'(^  V;   )       '>■  k^nS/  =  0 
XX         yy         '     /^     s   z 

1j  _  _    -  —     and  n       =  cm. 
From  formulas    (U.S)   ue  get   expressions   for  2  and  H  in  terms   of   the  new  V  and  './; 

\,;jir:i7.?..ctiely , 


^\  "  ^xz     *  ^^  ^  'V  ^^   "   •*'-"     -   ^^-■-^^' 


S3^  =  V         -  ik£    ',; 
•^  2        yz  X 


/     i       xz  r    y 


yuH^  =  W^^  .  ik/4  V^  (5.3) 


■^:ifc: 


*'    3  XX  yy  r    3  XX  yy 

The  continuity  conditions  are  the   same  as  before:      V.   \L   ■?- 7      .-   and     —  \I    must  be 

continuous   every\/here. 

Our  aim   is  to  express    tiic  functions  V(x,y,z)  and  \i{x,y,z)  by  a  complex 
integral   and  also  relate  V  and  M  to  the  source  \)f  the  electromagnetic  field,     've 
note  first  tliat   special   solutions   of  equations    (5.2)  are  obtained  by 

V(x.y.z)  =  v(z.c(,p)   e^^^^°^^//S^^ 

-,/  \        I  s     i^(=^  X  +£,y) 

'.'(x,y,z)  =  \r{z,<x  .p,)   e  *^ 

v/ith  arbitrary  parameters  c*  and  ^  if  v(z,Q*>,  A)  and  if  w(z,oi.    A  )   satisfy  the 
ordinary  differential   equations  ; 

^  (5.5) 

/X(^)      +  k^^^  v;=  0 

"^2  2  2  2 

v/here    0       =  £ i-<  ~  O^      -A      =  lu~  P     and  -orimes   denote  differentiation  v/ith  respect 
to   z.  '  '  '       '  * 
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All   the  follo\/ing  considerations  are  valid  also  in  case  of  conductinfl;  media(O"-^0). 
function  (.  has   onlv  tn  "h^   tp-hI  ar^or^   Y.\r  +Vio  r.^r-,v^i.,^  ■p„_„4.j -_  c*_  £       4  rfcr   ^ 


The  function  £  has   only  to  be  replaced  by  the  complex  function £"=  ^-  Vt^  i. 


uT 
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If  we  no\/  think  of  —  ?.r.d  -h—  as   the   direction  cosines   of  a  ray  with 

n      n 


respect  to  the  x-  and  y-axes  .tl-ien,  v/ith  so::ie  justification  ^-'e  nay  call  these  special 
solutions  plane  \/aves.   Indeed  "by  solving  the  second-order  equations  (5,5)  in  any 
region  \;here  S   and  i<  are  constant  the  solutions  ha.ve  the  form 

^ik(  ex  X  -y   /3y  -  /i^  -0^2  -/3^  2) 

which,    foro<      +p       =   n   ,   are  ordinary  progressing  plane  v/aves.      However,   v/e   do  not 
exclude   the  case<^     -*■  y      >  n   .    in  the   follouing.and  refer  to  these   solutions  also 
as  plane  uaves. 

Because   eqxiations   (5.2)   are  linear,    more  general  solutions  can  "be  found 
■fay  superposition,   naxaely,   by   integrals   of  the  type 

V(x,y.z)  =    Jj^{z',  &.,B)   e^^'^^"  ^  ^P^^   dc<  dp  . 

^^  '  .        (5.6) 

V(x.y..)  =     ijv;(2;c.,^)  e^^(^^^/5^^  dcdp. 

assuming  that  differentiation  under  the  integral  sign  is  permissible  for  the 
functions  v  and  v;,  xrith   \iiich  v;e  shall  deal,   Eq'uations  (5-d)  sho^;  the  forn  in  which 
ue  desire  to  express  the  solution  of  our  propagation  problem.   Hcv/ever,  they  do  not 
represent  the  final  an3\.'er  to  the  propagation  problem  for  \/e  have  yet  to  relate  V 
and  'd   to  the  source  of  radiation. 

li'e  nov/  assume  that  at  t^ie  point  x  =  y=z  =  0,'a  source  of  radiation 
is  located.   Some  kno\/ledge  of  this  source  i.rast,  of  course^  be  available,   w'o  suppose 
that  the  source  is  given  as  a  superposition  of  plane  \;aves  in  a  homogeneous  medium 
I     ,  u    ■.     More  definitely, ue  suppose  the  soxorce  given  by  the  complex  integrals  -^ 

z>0:    V=  j|?^(c<,/3)  a^^^"^^"^^"  ^V  °  '^'   -^  do<  d^ 


(5.7) 


(2  _  2    s 

ik  (  «  X  -t-  ,8  y  -  z  v/  '   ^ 


z<0:    W=  jj  G^Coc./J)  e^^^^^^"".-^-^^-  ^V'"         do(  d/^ 


-^'The   discussion  of  integral  representations   of  v;ave  functions    on  pp,   36I-36U  of 
Stratton,   J.  .a.,:      Slectromagnetic  Theory,   may  be  helpful  here. 
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in  which  n^  =  J'i^^^,    p     =o<     •*■  ^   ,  k  =  .^  ,  and  the  1^,   F^,  G^  end   G^  repre- 
sent the  amplitudes  \;ith  \;hich  the  plane  v/aves  enter  into  the  combination.   These 


integrals  converge  if  \/e  make  the  arjreement  that  the  sign  of  "2"  =  (/n.  -p        is  al\mys 
to  "be  chosen  such  that  the  imaginary  part  ol  "  is  positive.   If  o   is  real,  the  + 
sign  is  to  be  chosen.   The  ranges  of  integration  in  formulas  (5.7)  are  from  -  co  to 
00  foro;  and  i^  .  Generally  t\/o  expressions  are  needed  to  represent  a  source  "because 
the  source  itself  has  to  "be  represented  mathematically  "by  a  singularity. 

As  an  important  special  case  \ie   mention  the  dipole  oscillating  in 
direction  of  the  z-axis.   In  this  case 

F,  =  ?^  =  — 


^   2   ^  (5.S) 

0       ^ 


G,  =  G^  =  0 


The  complex  integrals   become 


iKn  r 
o 

V  =  ^  .,...       :      V  =  0  (5.9) 

J.  K   1 


The  proof  of  the  result    (5.9)  use-s  a  basic   result    in  propiagation  theory. 
If  v;e  v/rite  the   first   and  second  expressions  of   (5»7)  wi'th  the  values   of  F^    and  F_ 

given  by   (5-S)  \ie  obtain 


/2  2^2, 


\-Ie   change  from  c<, ,  A  coordinates  to  polar  ccordir^tes^  /J,  (p  ,   Since 
cn  and  Q   range  over  the  entire  plane  \/e  can  likeuise  cover  the  entire  plane  by  allow- 
ing i^  to  range  from  0  to  2n  vmile  P   ranges  from  0  to  oo^  Under  the  transformation 

^  =  p  cos^p* 
/2  =  p  sin^ji; 
with  the  chrnge  in  the  element  of  area  from  dtx  d/i  to />  d^  d<^  ,  expression  (5. 10) 
becomes 

V=    I       ^ ^ik(  ^;cosV5  X  4-^sin(^  y  +  zv/no  -/O   L^^^,, 

2iT  ^n     -  (>^ 

0   ' 

V/e  now  let   £>  =   n        sin  Z/     ; 
then       d/i  =  n     coz  p  d.'S/ 
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The  range  of  integration,  for  "^  must  nou  "be  nade  complex  to  cover  the  range  from 

0  to  00  for  O  ,  A.S   P  goes  froLi  0  to  n  ,  sin  y/ goes  from  0  to  —  ;  as  «>  goes  from 
n  to  00  ,  "^  must  range  from  p  to  ^  -*■  i  co  .   V/ith  this  change  of  variable  the  ex- 
pression for  V  ■becomes 


2"   ^°°  n 
0 

0 


ikn  (x  cos;/Csin'yy  +  y  sin.'i/sin  p+  z   cocy/  )   .  .,,  ,,,,  ,  ,„ 
e   0      1^         r        •>  f  r  1      '   3in  y  d^  d^ 


Excet)t  for  minor  differences  in  notation  and  conventions  as  to  direc- 

lU 
tions  for  measuring  angles  this  e:>qoression  is  the  sane  as  the  one  given  in  Stratton 

iknor    \\\ 

for   r^ . 

.  ikr 

The 'question  arises  as  to  hot/  ■i/ith  a  given  antenna  system  one  could 
knov/  the  functions  V  and  W  of  such  a.  source.   Theoretically  the  process  could  be 
this:   one  might  l-oio'J  the  eluctromagnetic  field  given  off  by  the  source;   in  par- 
ticular, one  might  knov;  S^  and  H   .  3y  using  equations  (U.  6)  functions  V  and  W 
can  be  obtained  uhich  are  not  unique  but  \/hose  arbitrariness  v;ill  not  affect  the 

ultimate  solution.   Given  V  and  V/,one  can  invert  the  double  Courier  transforms  in 

15 
expressions  (5.7)  ^o  ^i^*c^  ^]_.  ^2*  "^1  ^'^'^   ^2 

?or  some  antenna  systems  one  kno'..'s  the  3,  and  H^.   For  example,  the 
field  of  an  array  of  dipoles  is  kno'./n.   Likewise,  field  expressions  for  a  paraboloid 
fed  by  a  dipole  have  been  obtained  theoretically  in' the  past  v/ith  fairly  good  agree- 
ment between  theory  and  practice.   Tor  antenna  systems  such  as  horns  the  field  ex- 
pressions  are  also  kno\/n  approximately^  ?or  paraboloidal  antenna  systems  v;ith  more 
complicated  feeds  some  theoretical  \;ork  has  been  done  v/hich  gives  some  knov/ledge  of 
the  field.  ?or  example,  Radiation  Laboratory  Seport  T-9  and  7^2-1  by  R.  C.  Spencer, 
relates  the  illiunination  supplied  to  the  aperture  of  a  paraboloid  by  the  feed  \;ith 
the  field  pattern  produced  by  the  entire  assembly.   It  is  also  possible,  in  some 
cases  v/here  the  radiation  patterns  of  anterjia  systems  are  knovm.  only  by  measurement, 
to  obtain  the  functions  V  and  '</.  Measurement  usus.lly  produces  the  square  of  the 
field  intensity.   In  those  cases  v;here  the  polarization  is  vertical,  say,  the  square 
root  of  the  measured  value  is  3  .  Ho\/ever,  since  E^  would  be  known  only  riumerically 
it  would  be  necessary  either  to  fit  a  function  to  the  data  and  then  proceed  as  above 

■*■  'Stratton,  loc.  cit.,  p.  57S,  formula  (27). 

Compare  formulas  (12)  and  (I3)  on  p.  3^3  0^  Stratton,  loc.  cit. 
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or  one  might   find  V  and  ':J  tlurough  equations    (U.6)  "by  sor.e  n'omerical  process   and.  then 
invert   fornulas    (5«7)    'to   obtain  ?,  ,   F„,   G,  ,    and  G   . 

Ve  return  nou  to   one  najor  protlem  which  is   to  determine  the  electro- 
magnetic  field  radiated  "by  a  general  source   (5»7)   i^  i^    is  located  atx  =  y=z  =  0 
in  a  medium  of  stratified  character  given  "by  c  (z),    A-(z);      that   is,    \/e  laist   some- 
how involve  the  source   in  the  equations    i^.o)   for  V  and  U, 

6.   Propagation  of  an  Arbitrary  Source  Located  Belov;  a  Half-Space, 

"if.-y'M  As  a  step  in  the   direction  of  a  :.iore  general  propagation  problem  and  as 

an  indication  of  the  method  of  relating  a  source   of  radiation  to  the  field  propa- 
gated let  us  consider   the  problem  of  a  source   of  radiation  located  below  a  half-space 
of  arbitrary  electromagnetic  character. 

Specifically,   a  source  of  radiation  is  located  at   the  point    (0,0,0).     A 
stratified,   but   otherwise  arbitrary,   medium  is  supposed  to   exist   in  the  half-space 

"z  >  0  and  a  homogeneous  medium  is   supposed  to   exist   in  the  region  z   <  0.     The  source 
is  known  to  us   in  the    form  of  equations    (5.7)'      "'e  shall  consider   in  the   follov:ing 
pages   only  the   function  V  since  the   corresponding  considerations   for  \I  only  involves 
replac  ing  £,  by  yU  « 

We   shall   approach  this  problem  by  using   the  fact   that    the  "source  can  be 
built  up   of  plane  v;aves  and  begin  a  treatment   of  the  above  problem  v;ith  the   source, 
for  the  moment,    a  plane  wave. 

'tie  assume  that  S   and.- /<  have  constant  values  £    ,     Aj     for   z   <  0,    and 
that  £-^£,-,,     yU— >yuii^z  -^  oo  .     A  plane  \/ave 


V(x,y,z.t)  =   e  '^  .  e   ""*  (6.10 

is  assumed  to  be  incident  from  below.   The  electromagnetic  field  established  in  the 
region  z  >  0  under  the  influence  of  this  incident  plane  v;ave  then  can  be  obtained 
from  a  function 

V(x,y,z,t)  =  v(z,o{  ,  ,^5  )  e        r*"e  (5.U) 

V'here  v(z,a,^)   satisfies   the   ordinary   differential   equation 

iC-^)    +  k^-j'^^  =  0  (5.5) 

v/here 


•^   =v/£/U  ^os^  -f       =.     /n2  _    /?■ 
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'•  'j'e  are  led  to  seek  a  solution  such  as    (5.^-l-)   ir^  the   region  z   >  0  because, 

physically,  the   effect   of  the  non-horuogeneous  nediuia  should  te  on  the  liehavior  of  the 
v;ave  in  the  z-direction.     Put  nathenatically.v/e  raay  say   that   if  v;e  can  find  a 
V(x,y,z,t)  as  given  ty  equation  (5-^)  v/here  v(z,c<ji3)  satisfies   equation  (5.5)'.-'^'e 
can  obtain  an  eloctroaagnetic   field  \/hich   satisfies  Haxv/ell's   equations   and  \/hich 
can  "be   interpreted  as   the  alteration  of  the  plane  wave  arising  frosi  the  non-hoao-    . 
geneous   layer, 

V^e  assume,    at   first,   more   specifically,   that  £•  =  i       U  =  Lk     for   z   >  a, 

the  quantity  a  being  chosen  arbitrarily  large.  q     ^  i  ) /"I 

In  the  hor.-iogeneous  parts   z   <  0  and  z>J(         —   —   _^    —    — 
of  our  niediun  v(z)  =  v(z,c«.  ,^  )  r.ust  have  the  fora  ^^  ~'        p] 


ik  '2^   z  -ik  "iT  z  T" 

z<0{v=e°+Re  ° 


>/  :  ^  =  T  e 


-;>X 


^.. 


'/"' 


That  is,  in  the  region  z  <  0  \/e  allow  for  the  presence  of 
an  incident  wave  of  unit  amplitude  and  a  reflected  wave  of  arrplitude  R,  while  in 
the  medium  z  >/ mq   allov;  for  a  transaitted  wave  of  aciplitude  1,  The  anplitudos  R 
and  T  deterLiine  the  reflectance  and  transaittance  ;of  the  inhonogeneous  nedium  v/ith 
regard  to  the  incident  plane  wave  (6,1),   The  values  of  a  for  z  .<  0  and  z  >/<  are 
r  and  r, .respectively, 

'.e  refomulate  the  boundary  conditions  (6.2)  by  eliminating  R  and  T 
with  the  aid  of  the  derivatives  v'(o)  and  v^{X   ), 

We  find  v(o)  and  v'(o)  from  the  first  of  equations  (6.2)  and  eliminate 
R.   We  then  find  v ' ( x  )  and  v(x  )  from  the  second  of  equations  (6.2)  and  eliminate 
T.   This  gives  ^ 

v»(o)  +  ik:2r  v(o)  =  2iky 

0  0 

(6.3) 
v'(i')  -  ik/^  v(i)  =  0 

The  plane  wave  problem  now  is  to  find  a  solution  of  equation  (5.5)  which  satisfies 
these  boundary  conditions. 

We  can  reduce  this  problem  to  the  problem  of  finding  a  solution  of  a 
certain  differential  equation  of  first  order  of  Riccati's  type. 


on 
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Me   introduce  the  dinensionless  function 

1  16 

Since  j-  v'(z)  and  v(z)  are  continuous  functions   we  conclude  that  ©(z)  nust  te 

continuous  for  all  z  even  if  £(z)  and  //.(z)  sho'old  liave  finite  discontinuities. 

We  introduce  -r—  =  ik  v  0  in  equation  (5-5)  ^'■^   obtain  an  equation  for  ©, 

c 
nanely  ,  P 

©  =  ik(4--6e^)  (6.5) 

i.e.  a  siuple  Hiccati  equation.   The  second  condition  (6.3),  together  \/ith  equati 
(6.^)  gives  us  the  "boundary  condition  for  ©(z)  at  z  =  X  . 

Letting  J(  —>   co  we  recognize  'S(z)  as  the  solution  of  the  Riccati  equation  (6.5) 
\rhich  satisfies  the  condition 

Qiz)    ->  — ^  if  z  — >  00  (6.6) 

^1 

The  first  condition  (6.3)  together  \;ith  equation  (6.U)  can  nov/  te  used 

to  express  v(o)  in  terns  of  0(o).  Also  by  the  first  equation  (6.2),  v(o)  =  1  +  R. 

These  two  facts  together  give 

2  V 
v(o)=l+R=    ^      ^°, .  (6-7) 

where  ©  =  0(o)  but  <£   is  the  value  of  £  ii^  the  medium  z  <  0  \/hether  or  not  £  (z) 
is  discontinuous  at  z  =  0.  The  reflectance ,R  itself,  is 

^0 

a=    °  , (6.S) 

^o 


^0 
Finally  the  solution  ,v(z)  itself,  is  found  by  quadratures  from  equation, 
(6. U), namely 

log  v(z)  =  y   ik£e(z)  dz  +  C 

ik/^<£G(z)  <i2 
v(2)  =  D  e  ° 

^6.  If  ;;e  select  a  continuous  solution  of  v  for  0  5  z  i  )i  of  equation  (5-5)»  which 
meets  the  boundary  conditions  (6.3)  then  v(z)  v/ill  be  continilous  every\/here.   Then, 

s  V  /",  2  ^2 

by  equation  C5.5). -7-  is  expressible  as  /  JL__^__J1  dz.  Even  if  £  has  finite  dis- 

continuities  this  integral  is  a  continuous  function  of  z  for  all  z. 
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Evaluating  L  at   z  =  0  gives 


v(z)=v(o)      e 


z 


ik   /£©(z) 

=  vCo)      e 
7 
or,  using  equation   (6.Sj, 


2   <■  ik  y    ci  ©  ( z )  dz 

■*  0  0  ^-^0 

Returning  now  to  equation  (5»^)  ^^d.  onitting  the  time  dependence  ve  obtain  the 
function 

.  /.z 

2    ^''  ik(cKx+/^yW    £   (z)e(z)     ds) 

V(x.y.z)=     -——^-^     e     ^-  ■  °  -^  (5.10) 

»  0       c  0   C^  0 

This  function  is   the   solution  for  V  of  our  plane  v/ave  proolem  in  the  non-honogeneous 

nedixiQ   z    >  0.      3y   sulDstituting  the   first   of   equations    (6.2)   in   equation   (5.^)f   v/ith 

the  value   of  R  given  "by  equation  (6.8),  v;e  ootain  the   solution  for  V  in  the  region 

z   <  0. 

Just  to  "be  clear  on  the   order  of  events,  were  an  actual  calculation 

to  te  undertaken,    one  v/ould  "begin  "by  solving  equation  (6,5)   foi'  G  (-)   subject   to 

the  boundary  condition  (6.6).      One  then  obtains  ''/(z)  fron  equation  (6»9)   and  then 

finally   equation   (6,10).      The  solution  in  the   region  z   <  0  is  then  iraiiediate  fron 

equation   (6.8)   and  the   first   of  equations    (6,2).      The  theory  also  covers   the   case 

v;here  k    is  finite,    in  v;hich  case  three   functions,    one  each  in  the  regions   z  _  0, 

Oizi    \  ^   z2    X,    represent  the   complete  solution  for  V.      In  this  case  one   needs 

also  the   second  of  equations    (6,2)  \/ith  T  =  v(  O.     Also  the    solution  (2>(z)   of 

equation   (6.5)  auot   satisfy  the  boundai-y  condition '3  (  X  )  =     — — -  , 

^1 
ik^^z 
Since,   by  equation   (6, 1),  v(z, c< ,  5 )   =   e  in  the  case  of  a  plane 

v/ave  passing  through  a  homogeneous  ;nediuja,   it  follov/s   fron  the   definition  (6.U)   of 

0(z)   that,  in  an  unbounded  homogeneous  raediun,0(z)  =  __2.     every\/here.      If  v/e 


now  compare  equations    (6. 10)  with   (b.l),    to  which  equation  (6. 10)  reduces  when  the 
aediuLi  is   everyvmere  homogeneous,   we  see  tlie    significance  of  the  function  Q>{z), 
Through  the  expression 

cx.x+/3y-^/     £(9(z)     d? 


2U. 


it  determines  the  distortion  of  the  wave  fronts  c^  x  + /3  y  "^  oz  due  to  the  non- 
hoQOgeneous  mediua  and  through  the  expression 


=■0   '-o 


it  determines  the  change  of  anplitudc  of  the  undisturbed  wave. 

A  similar  result  is  obtained  for  the  function  v;(x,y,z); 


W(x,y,z)  = 


2/. 


ik  1  ex  X  -»-  /^  y  +/  U-Q   dz 


[6  y 


0  / 


where  6(z)  is  the  solution  of  the  Hiccati  equation 

2 


(6.11) 


0'  =  ik 
which  satisfies  the  boundary  condition 


^ 
4 


0(z)   -» 


/^G' 


(6.12) 


A'l 


if  z 


00 


->A 


£ 


o^/^i 


6.1  The  above  results  allo\;  us  to  give  immedi- 
ately the  formal  solution  of  the  problem  of  propagation 
\/ith  which  v;e  began  article  6,  namely,  the  propagation  of 
an  arbitrai'y  source  located  belou  a  half-space.   The  source 
of  radiation  is  placed  at  the  point  x=y=s=0,   A 
Gtiatified  medium  is  c^oncidered  './hich  is  homogeneous  for 
z  <  0.   The  radiation  pattern  of  the  source  in  a  homogeneous 
medium  uith  characteristics 
integrals  of  the  form 


^0- 


La.     shall  be  given  by 


;■    r  ik(o<.x  +/5y  +  Yz) 

V  =     .'    /    r,(o<,/6)   e  °       do(d/3 

J  J 


for  the  region  above   the  source. 


ik(o<x  ^^y  +  V^z) 


(5.7) 


dc< 


id.fB 
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If  this  source  is  placed  teloi.'  the  stratified  mediurn  \/hich  "begins  at 
z  =  0,  then  the  associrted  electromagnetic  field  in  the  region  z  >  0  is  dctcrnined 
"by  the  tv;o  functions 

2 


VCx.y.z)  =  2 


V/(x,y,z)  =  2i 


^0 


/r 


ik  fcK  X  +5  y  +/  £  ©  dz  1 
-3^?,(o..^)e    L     r     0      J  ^^^ 


^  ik  |o<x  +/3y  V  /<o  dzl 


(6.13) 


"here    y^,     £q.   and   ^^ 


^0  ■*■  A 

are  the  values  in  the  homogeneous  region,  and  t'here  0(z) 
and  0(z)  satisfy  the  Riccati  equations 

©•(z)  - 


2ni  ,  2_ 

A      C 


LO^) 


o'(z)  = 


X 


'N  J   -  ;.o^, 


(6.1U) 


and  "boundary  conditions 


©(O  ^ 


r. 


^1 

G(z)  -»     — i 


>  if 


•^     00 


The  ranges  of  integration  foro<  a.nd  A  in  formulas  (6.13)  are  understood  to  "be  the 
same  as  they  are  for  the  source  in  formulas  (5«7)'   --e  components  of  the  electro- 
magnetic field  itself  are  nov;  o'btained  "by  equations  (5'3)' 

The  Justification  of  this  general  result  is  simple.   The  exiDressions 
(6.10)  and  (6.11)  \.'hich  v/e  have  ohtained  satisfy  the  requisite  differential  equation 
(5.2)  "because  equation  (6.3)  aiid.  its  analogue  for  the  magnetic  field  satisfy  equa- 
tions (5«5)«   Bxpressions  (6.10)  and  (6. 11)  then  are  the  solution  of  the  propagation 
problem  of  this  article  for  a  plane  \;ave  source.     If  v/e  multiply  these  expressions 
for  "V  and  U  "by  constant  factors  ?  (■^,|^)  and  G^(c\',|3  ),  respectively,  for  fixed 

o^  and  ^Q  ,  we  still  o'otain  solutions  of  the  "basic  equations  (5-2)  for  V  and  \l.  Be- 
cause equations  (5.2)  are  linear,  any  sum  of  solutions  finite  or  infinite  (provided 
convergence  criteria  are  satisfied)  is  also  a  solution.  If,  however,  the  source  is 
no  more  than  a  superposition  of  plane  \/avcs,  the  propagated  field  of  this  radiation 
source  in  the  non-homogeneous  medium  is  no  more  than  the  sum  of  the  solutions  for 


26. 

Uve  indiv-i&ial-pla.n6  v/avec-, — i^Hs^Aa-s-4-Sr-I3-)~s-'^^''^^:>"~~^<^^  1  ou 

ro^  the  individual  plane  waves.  Fonnulas  (6.13)  amount  to  no  more  than  the  summa- 
tion of  an  infinite  number  of  plane  \/aveG  each  \;ith  itc  anDroTDriate  amplitude  factor 
F.Co^,  (2)  (or  G,(>>',/3  )),  and  hence  must  "be  the  solution  to  the  problem  of  propagating 

the  arbitrary  source  (5»7)  through  the  non-homogeneous  half  space. 

V/e  note  that  formulas  (6.I3)  are  very  general  for  they  cover  any  varia- 
tion in  E,  andyW.  f or  z  >  0, including  finite  discontinuities.   The  evaluation  of  the 
integrals  may  he  a  problem  of  come  difficulty  depending  upon  the  functions  ?,  and  G^ 

and  on  t   and /<..  Formulas  (6.I3)  are  not,  of  course,  the  general  solution  of  the 
propagation  problem  above  a  flat  earth  .  This  problem  \/ill  be  treated  in  article  f , 

6.2.  We  shall  consider  one  illustrii-tion  of  the  theory  Just  developed.  As 
bef ore  ,an  arbitrary  source  given  by  erq^ressions  (5.7)  is  located  at  the  point 
(0,0,0).   The  non-homogeneous  medium  shall  consist  of  t\/o  layers  with  S  =  £.   and 

J^   =  ^     for  z  <  A,  and  <E  =  £.  andy^=  ^,  for  z  >  /  ,  ^  >  0.      Thus  our  medium  has 

a  finite  discontinuity  in  £  and  /-^ at  z  =  A  , 

In  accordance  with  the  general  result  of  article  6.1j\/e  must  first  find 
solutions  (3  (z)  and  ©(z)  satisfying  equations  (6.IU)  and  boundary  conditions  (6. 15). 
The  solutions  in  the  region  where  ^C=   V.  ,  c,=  £    and.  U=   /U,  „    v;hich  meet  the 

boundary  conditions  (6.I5)  (and  v/hich  therefore  must  be  unique)  are 


sr 


>i'   :   ©(z)  =     ^       .  0(z) 


^1 


-1  \  r-1 

In  the  region  z   <    X&{2.)  and  0(z)  must,    first  of  all,    satisfy  equations    (6.IU),  and 
be  continuous  with  the   solutions  for  z>Aatz=A.     To  solve    equations   (6.1U) 
when  Y' .  C.   ,   and yu.  are  constant  i/e  i/rite  then  in  the  form 

©'(z)  _       2ni 


v/hich  are   immediately   integrable  by 


e'(z) 

2ni    y2 

r 

k 

X     ^ 

, -1          _     /         dx 
tan  h         X  -    /  


o      1  -  X 
lo  meet   the  continuity  condition  we   require   that 

©(i?)=     ^     and     G(2)=     II 
^1  Al 


27. 


Carrying  out  the   integration  and  detercination  of  the  arbitrary   conctnnt  by  this 
condition  ,ue  obtain  for   z   <  JC  , 


0(2) 


1       .     ^0    .        2n  f  Q         V    ^ 
-  X  -T—  tan  -^  (/  -  z)  Y^ 

'-o  o 


r. 


a, 


1  -  i 


^1    4< 


,  e(z)  = 


^_i_:^tan|IIV'(/-   z) 


<^l        >C. 


tan 


?(J?-  z)  Y. 


X 


1  -  i 


<1   /' 


:( 


tan^rCi'-  z) 


./«1     To  ^       ° 


To  use  formula   (5,13)   ^'<i  need  also  0(o)   and  O(o),\;hich  axe  inmediately 
obtainable  from  these   last  fonailas. 

The  fornulae   (6. 13)   simplify  to 


z    <  ,^  :      V  =  V     +  V^ 
0  1 


■■uU.O'         z   >  X  :     V  =  V, 


w  =  :•!    +  \L 

0  1 


\i  -  \K 


\/here  V     and  \1     are  the  t;;iven  incident  v;aves    , 
0  0 


V     = 

0 


ik(cx  X  +/2y  +Y  z) 


w       — 
0 


Gi(<^,{3  ) 


ik(c^x  +.5y  +  £j^  z) 


do^d.5 


V,    and  ■<■     ,the  reflected  \;aveG  , 


-f ^    Fi(o<.^)   e        '  ^  °  dCd^ 


■•^1  = 


0 

/^o 


ik 


^. 


1 


oix 


^/3y  +   Y  (2/^-   z    )] 


dOdB 


/^o  /^l 


2£. 


i-v^us-i^jV-  and  :/    ,the   transmitted  \/aves  , 
!    I    d  e 


^1 


.(^.^) 


ik  [=<z  +-5y  ■»-  Y'3_(z  -  b] 


di?<d(3 


ra' 


?  -—2.  p       "  °  -  -I 

JU.^  ik  kx  -^.^y  ^    -iT,  (z  -.^)| 


ri 


In  case  x^C^C./S) 


'/^ 


2n7n^-os2-yO,2 


2n    V 


we  have  a  dipole,    oscillating  in  the  direction  of  the  £-axig. 


7,   Propagation   of   an  Ar"bitrary  Sourcs  Located  3et;;een  Tv;o   Ealf-Spaces« 

Oar  raajor  goal,    naiaely^    the  propagation  in  a  non-homogeneous  atmosphere 
of  ■ultra-high  frequency  radiation   from  a  source   located  alDove  the  earth  \fill   nov;  "be 
approached  through  a  more  general  proolem.      V.'e  shall   suppose  the  source  to  he 
located  at  the  origin   of  a  rectangular  coordinaoe   system  and  allov;  the  half-space 
above  and  belo\/  the  so\irce  to  have  ar^i^  contiroious  or  finitely  discontinuous  varia- 
tion oi  S  ,  LL  ,   and  O  .      './e  uay  then  specializevthe  medium  below  the  source   so  that 
it   represents  the  characteristics   of  the  earth. 

As   in  the   treatment  of  the  problem  of  article  6  ve   shall  first   allov;  the 

source   to  be  a  plane  uave  and  build  up  the  arbitrary  source  from  plane  trcives.      We 

consider  then  a  medium  \mich  is  homogeneous   in  a  small  strip  -  ^  <  z   <  S   but  other- 

v/ise   inhomogeneous.      Let  us   imagine  that   the  plane  z  =  0  is   the    source  of  tv;o  plane 

waves   symmetrically  oriented  to   the  plane  z  =  0. 
>2   A 


r 


In  a  homogeneo\is  medium  these  emitted  plane  waves  shall  be 


ik  ^f   z 


^1^ 


for  z  >  0 


^:Sm$7^M^->  / 


and 


-ik  ^  z 
Fg  e        for  z  <  0, 


■■"^ 


29. 

In  order  to  deteriaine  the  influence  of  the  medium  ue  introduce,  as  auxiliary  quanti- 
ties, two  reflection  coefficients  H,  and  R  .   R.,  is  the  reflection  coefficient  of 

the  upper  medium  z  >,5  if  the  lo\;er  medium  z  <  d  is  homogeneous,  R  is  the  reflec- 
tion coefficient  of  the  lo\.'er  medium  if  the  upper  medium  is  homogeneous. 

\/e  utilize  earlier  theory  \.'herein  the  solution  of  an  elect^o^.^agnetic 
field  for  a  plane  uave  through  a  medium  vpr/ing  in  the  z-direction  only  is  reduced 
to  the  determination  of  v,  as  in  equation  {5. 5),  and  make  the  substitution, as  before: 


0(z)  = 
v/hich  gives  as  "before 


1    v'(z) 


ik   £  vCs) 


e'  =  ^(  |-  -  £G^)  (6.5) 

The  physical  conditions  are,  of  course,  different  from  those  of  article 
6,  since  \/e  no\/  have  media  of  variable  characteristics  above  and  belov/  the  region 
-  S  ^  ^  '^  E     ^^^   '''^■^^  incident  v/aves  are  in  t\/o  directions.  Moreover,  t;ro  different 
solutions  of  equation  (6.5)  a-^«  needed  as  compered  \;ith  one  in  article  6,for  tv;o 
boundary  conditions  nov;  must  be  met,  one  at  2  =  +00  and  another  at  z  =  -  00  .and 
the  solution  of  a  first  order  equation  can  satisfy  onl;,-  one  boundary  condition.   Let 

QAz)   andO  (z)  be  t\/o  solutions  of  equation  (6.5)  which  shall  be  deter::iined  by 
the  conditions 


e,  (z)  ->  4^  if  z  ->  +  00 
^       ^1 

9o(z)  "  -  v-^  if  z  ->  -  oov 


(7.1) 


V/e  may  no\/  take  advantage  of  the  theory  applied  in  article  6  to  obtain  the  reflection 
coefficient  at  each  boundary  z  =  b  and  z  =  -a    ,     As  R.^  and  R  are  defined  here,  R^ 
can  be  determined  for  the  plane  wave  directed  up  and  R  for  the  one  directed -down, 
'jithout  regard  to  the  other. 
Hence  (cf.  formula  (6.8))  /■ 


^-©,(0) 

R^  = 

%  *  ®=(°> 

1°  *  e,  (0) 

^  0 

2 

e:  -  "^'°' 

Si^ 
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l?ov;  if  the  nedioim  "below   z  =o    '..'ere  hoiTiOgeneous,    the    field  for   z   <$ 

ik  '  ii^   z 
\YOuld  consist   of  the   incident   field     7^    e  and  the  reflected  field 

ik  iT   z 

R^    ?,    e  .      Likev/ice  if  the  nediim  atove   z  =  -  o  \'&re  homogeneous,    the   field 

-ik  1^     z  -ik   "sf   z 

arising  from  Fp  e  v;ould  consist   of  this   ter^n  and  Il„  ?     e  .      Houever, 

the  presence  of  "both  plane  ifaves   and  both  ooundf.ries  means   that    in  the  region   z   >  0 
v;e   shall  have   radiation  directed  toirards   the  upper    ocundary  which  is   eaittcd 
directly  frou  the   source  of  uagnitude     F,  ,plus  radiation  reflected  from  the  lov;er 

boundary   z  =  -S'.of  ■uehuo'Si  lia^rnitudo  A  ,   All  of  this  will  be  partially  reflected 
froa  z  =  +0    .      The  reflected  radiation,   of  aiiiplitude  3  say,   will  be  directed  do\;n. 

Hence  there  vn.ll  be   radiation  directed  towards   z  =  -  0  v;hich  consists   of  the  source 

-ik  Y   z 

Fp  e        plus  \;hat  is  reflected  fron  the  layer  z  =  S   of  magnitude  B.   Reflected 

from  z  =  -S  will  be  radiation  of  magnitude  ^n^'?  "*"  "^^  "hich  in  turn  must  be  A. 

Likewise  reflected  from  the  upper  boundary  z  =  -^  v/ill  be  radiation  of  magnitude 
R,  (F^  +  a)  v/hich  i.Tust  be  3. 

It  follov;s  that 

From  a  purely  mathemf-tical  standpoint  v;e  aay^  say  .that  in  the  homogeneous 

strip  -O  <  z  <  ^  the  solution  of  the  differential  equation 

\ 

£(:^)'.k2  ^2^  =  0 


can  take  the  form 


ik   X'    z  -ik    ^^  2 


v^  =   (F^  +  A)  e  "     +  3  e  ,      z   >  0 

,^^  ^o"  -ik   r  z 

v^  =  A  ■'+  (F^  +  5)    e  °    ,      z   <  0 


(7.U) 


v;here  v,    and  Vp  are  two   distinct    solutions*- 
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Then 


1  +IL 

(7o) 

1  *H 

V  (  0 )  =  ——   ( ?  +  ?  ?  ) 

""S^^    1  -  R^R^   2    11^ 

To  obtain  the  colutions  ir.  thj  riou-honogeaeous  uedia  ue  no'.;  treat 
V  (z)  and  v„(z)  as  in  article  6.   './e  repeat  the  steps  fonaerly  made  to  obtain  equa- 
tion (6.9)  from  equation  (6.U)  only  the  v^(o)  here  replaces  the  v(o)  equation  (6.7). 
'./e  repeat  the  process  once  r.ore,  this  tLv.t;  using  Vp(o)  to  replace  v(o).  ''e  then 
substitute  each  result  in  equation  (';>/4).  lIo\/  ci  is  arbitrarily  snail.  Hence  \7e 
Eay  say  that  .  •■ 


z 


1  -Ml^ 
<  0  :     V^Cx.y.z)  =     ^_^^       (?i  -  ^2?^) 


ik(o(.x  +  j^y  +  £02  dz) 

wo 


12 


(7.6) 


l-R^ 


ik(c<x  +  fiy  +j    £0^    ds) 
>C:       V^(x.y.z)=     "Tir^^   ^^2  "  Vl>   ^ 

This   is   the  solution  of  our  problen  of  propag^on  fot  plane  v/ave  sources*      If  the 
source   is   syixietrical ,    i.e.,    if  ?^    =  ?„  =  F  =  1  then  our  result    siciplifies   consider- 
ably by  introducing     O,(o)  and  0p(o)    in  (7.5)  by  aeans   of   (7-2).     '..'e  obtain 


^    <0   :   ^2Ury.0=      ^^(,)   ^   ^^(,) 


2'^o/£  i^(o;xV  ^ay  +y    J.  Q^   d2) 


(7.7) 


2/^,^  ik(c5'vX  +   py  -t- j      £G3_  dz) 

z   >0   :  V^(x,y.z)=     q^(,)  _^3^(,)       e  ° 

Thus  if  the  source  consists  of  two  plane  ua'^'es,  one  of  amplitude  F^ 

directed  upv/ards  and  one  of  amplitude  ?p  directed  do\mv;ard,  then  the  effect  of  the 
.-nediua  on  the  \/ave  is  given  by  formulas  (7.6)  ^:ith  corresponding  expressions  for  ':!, 
If  the  two  plane  uaves  have  ax-iplitude  1  then  the  effect  of  the  iijediun  is  given  by 
expressions  (7.7)  uith  corresponding- expressions  for  V/,  The  i'ull  field  expressions 
for  plane  waves  in  the  non-homogeneous  media  are  then  obtained  by  means  of  relations 
(50), 


32. 

V.'e   turn  nov;  to    the  major  proolen  of  a  general   source  of   radiation  placed 
at  X  =  y  =   z   =  0  v/ith  the   media  atove   and  belov;  non-homogeneous   to   the   extent   of 
being  stratified.      ?or   this  more  general   source  v;e  reason  again  as   in  article   6, 
The  character  of  the   source  must   ce  laicwn  to  us  either  through  the    functions 
F(c<, /S)   and  G(c<,  (3)   or  through  Y  and  \I  for  a  homogeneous  medium  from  \;hich  \/e  may 
find  P  and  G  through  relations   (5.7)  "by  inverting  the   Fourier  transforms. 

Again  we  argue  that   since  our  general   souree   is  a  superposition  of  plane 
vaves  each  modified  by  the   amplitude  factor  r,-.       .),    the  suTjerposition  of  the  plane 
\;aves  modified  by  the  presence  of  the  non-homogeneous  medium  and  each  multiplied  by 
its  a^iplitude  factor  FCoC^j!^)  must  be  a  solution  of  Haxwell's   equations,    for  these 
equations,    or  the  equivalent   formulation   (5.2)  and   (5.3)   of  this  paper, are  linear. 

Specifically,    let  us  assume   that   the  source   is   symmetrical   so  that 
?j_  =  ?2  "=  F(o<,/o)   and  &]_  =  &2  "^  G-(c<. » ^  )    in  formula  (5.7).     Then  we   obtain  the 
functions  V(x,y,z)   and  v;(x,y,z)  by  superposir^^  '^lane"  waves  of  the    form  (7.7).   '^'^^ 
obtain 

rr  ik(cix +|3y +/  £  ej_dz) 

rr  ik(o<.x +py -j^/xe^dz)         ^'''^^ 

,,^^  rr  ik(o<  X  +  ,S  y  +/£  a^dz ) 

z   <0i^ 


v='.='-'>=Tf^/iiii^(7T^'^'P'^       :.     ^«^f.3, 


v/here  "^^(0  •♦•),  i^(o  +),    and  £(0  -*■)  are   the  values  approached  by  V,i/<.,    and  £  as   z 
approaches   0   from  above,    and  Y(o  -) ,  L^io  -),   and  £(0  -)  soce  the  values  of  ^f ,   u., 
and   £  as   z   approaches   0   from  belov/;    and  v;here   Q,(z),0   (z)   and  6^(z),   9p(z)   are 


the  solutions   of  the  Riccati  equations 


x^2 


0-  =    ~  i~ ^0^)  <7.io) 

A  <5 
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which  satisfy  the  "boimdriry  conditionG 


if    2    ->  +   00  (7.11) 


Q^{z)  -^ 


and 


if    s   — >  -   CD 


li 
/^i 

It    is   understood  that    here, too, the   rariges   of   integration  for  i?\  and  ^  are   determined 
"by  the   source. 

The  results   enoodied  in  forinulas    (7.S)   and  (7.9)   taken  in  conjunction 
v/ith  foroulas    (5.3)    constitute  the  najor  general  conclusion  of  this  paper,      "hey 
solve   the  prolilem  of  propagation  in  stratified  non-homogeneous  aedia  in  the   follov/- 
ing  sense.      One    stairts  with  the   nature   of   the   source,    that    is,    one   either  lciio\/s 
the   functions   F(c;',  ,  A  )   and  G(<^  , /3  )    or,    alternatively,    the   functions  V  and  '/  of 
relations    (5*7)   from  uhich  ?  and  C-  can  "be  obtained.     Physically  this  means    that   one 
knov;s  how  the   source   can  be    cuilt  up  from  jDlane  v;aves.     Granted  this,    formulas    (7.8) 
aiid   (7-9)   are  used  to   calciilate  V   and  \!  of  the  modified  field  and  then  the    field 
components  can  be   found  by  means   of   formulas   (5.3)* 

The  nature  of  the  inhomogeneity  oi'^  the  medium  above  and  belo\;  the   source 
is, of   course,    implicit   in  the  functions   O-,,    ^^  ,    Op  and  6     v.'hich  are   determined 

first  by  solutions   of  eouations   (7.10)  v/ith  the  boundary  conditions    (7. 11).      The 
character  of  the  inhomogeneity  is   restricted  only  to  the  extent   that  £  ,   Zc,    and 
CT  must  be  functions    (continuous   or  v/ith  a  finite  number  of  finite  discontinuities) 
of  z.    It  will  be  recalled  also  that    if  the  medium  has  conductivity,    £,  mist  be  re- 
placed by£-     — — - —  .      Solution  of  equations   (7. 10),   \/hich  are   of  the  Riccati  type, 
may  have  to  be  approximated.     The   ease   or  difficulty  v/ith  which  one  obtains   solu- 
tions  depends  upon  the  nature  of  the  functions  £•  and  JU. 

It  will  be  observed  that  the  result   for  V  and  'iJ  in   (7.S)   and  (7.9)  are 
complex  integrals.      The  evaluation  of  these  integrals   is,    in  general   not    simple, 
even  though  exact   or  approximate   expressions   are  :3iov/n  for    3   ,  Op,   6-,    and  0        for 

in  the   integrals  the  nature  of  the  F(c\,  rt)   and  G(o<,j6)  also  must  be  considered. 


Articles      7-^      to  7-^    '.rill   take  up   cpecial   cases   of  this   general   result 
These   special   cases  not  only  will  give   sozlc   indicp.tion  of  the    diff iciilties  net  in 
evaluation  "but   \;ill    come  closer   to  ir.ore  practical    siouat  ions. 

7.1   The   Dipole   Source  "betveen  l^on--F.G:nogeneous   St.ratified  hedia. 

The   radiation  of  a  dipole   oscillating   in  the   direction  of   the   z-axis, 

is   obtained  "by  letting  F(^,  f^'>  =     2^~V ^^^  G(oc,  iS)  =  0.      (See  eouation  (5.8) 

above.)      Forimilas    (7.S)   and  (7-9)    thus  yield  f or,  z    >  0, 

-  r     ik(oCx  +/3y  -^J^    £0^   dz) 

\-    -2FfTwjJ'    G,(o)    -G,(o) '^'/^         ^^'''^ 

W^  =  0. 


and  for  z   <  0 , 


p.     ik((Xx  +/6y  +/^  £©2  dz) 

As    indicated   in  tne   discussion  following  equation   (6.10),    for   the   un- 


bounded homogeneous  raediun  ,0,    =     ,  (-} 


00       ^      _  ^/ 


0 


2  L        ' 


^0  ^0 

V 

and  v/e   find  .,  \ 

ikn  r 
o 

ikr 
as    it   should  be   {See   the  proof   of   foroula   (5«9)/ 

7.2  A  General  Source   Above  a  Perfectly   Conducting  Ground  Surface, 


The  general   theory  of  article   ]   allov/s   fgr  any  type   of   (stratified) 
30US  medium  above   and  belO'.r   the   arbitrary  sou] 
specialize  this   to   the   case  uhere   there   is   a  homogeneous 


non-homogeneous  medium  above   and  belO'.r   the   arbitrary  source.      V/e   no\.'  seek  to 

,rl ^-A 


o 


/V<3/^^> 


layer  of  height  h  belov;  the  source  and  beloi.'  this  a  per- 
fect reflector,  v/hich,  in  practical  problems,  v;ould  be 
the  earth.   The  source  itself  is  assumed  to  be  at  the      ^w— ^^    ^  o/"//? 
origin  (0,0,0).  K^  ,yx/S 


35. 

To  obtain  solutions  for  V  and  ':!   for  this  problea  \ie   shall  need  ex- 
pressions in  the  three  different  regions  z  >  0,  -h  <  z  <  0,  and  z  <  — h.  The  third 
case  may  "be  neglected, for  no  field  exists  in  a  pe|*fect  conductor.   The  region 
z  >  0  is  still  unspecialized  teyond  the  fact  that  it  is  stratified.   We  can  "be 
specific  \;ithin  the  homogeneous  layer,   '..'e  roust  first  solve  equations  (7. 10)  for 
the  region  -h  <  z  <  0  subject  to  the  follov/ing  conditions.  Within  this  region 
^  .   IJ  .   and  0   are  constant  and  denoted  by  S  ii    ,i^tA   o     . 

•    A^  •  0/00 

V/e  note   that  within  the  perfect  conductor  the  dielectric  constant 


whatever  the  value  of  C\  and  A,    our  boundary  conditions,   uhich     replace  conditions 

(7.11).    are 

GjC-h)   =  0 

OgC-h)  =  CO 

The  solutions  If  equation   (7. 10)  \/hich  hold  in  -h   <  z  <  0  and  satisfy 
these   conditions  are: 


^Z" 


.    ^  0     ^  2n  -.y^ 

2-'  1^     tan     ^ 


GoCz)  =     i4r^     tan     41  /(z  +  h) 


0„(z)  =  1     — —     cot     ^-   <;•    (z  +  h) 

y^o 


(7.1^) 


In  particular 


^  0     ^         2n  ./■ 


02^0)   =  i     f-     tan     -^,h 
C-o         ^ 

These  four  quantities  ©  (z),  Q Az) ,    QAo)   and  0^(0)  must  now  be  sub- 
stituted in  formulas  (7.S)  and  (7.9)  to  obtain  the  field  above  and  belov/  the  soiirce. 
One  cannot  proceed  to  evaluate  the  double  integrals  because,  thus  far  the  character 
of  the  medium  above  the  source  has  been  left  quite  general;  that  is,  the  four  quan- 
tities Q-iCz),  Qi(z),  G->(o),  and  0  (0)  are  undetermined.  (Particular  media  above 
the  dipole  will  be  considered  in  later  articles. 

The  theory  of  this  article   cells    for  a  homogeneous  layer  belo),/  the 
source.   The  parameters  c  ,/^,  and  (J  in  this  layer  can  be  any  constants  followed  by 
an  arbitrary  variation  in  these  parameters  in  the  non-homogeneous  medium  above  the 
source.   The  restriction  to  a  homogeneous  layer  belo\;  the  soxirce  is  merely  for  pur- 
poses  of  illustration.  The  general  theory  of  article  7  pe^^^ts  any  arbitrary  medium 
'•  The  method  of  solution  is  the  sane  as  in  article  6.2. 
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"between  the  Gource  and  the  perfectly  conducting  earth.   If,  as  in  the  present  article, 
the  layer  "between  the  source  and  the  perfect  earth  extends  from  -h  _  z  £  0,  the 
"boundory  conditions  O  (-h)  and  Q„(-h)  =  oo  would  "be  retained, 'cut  equations  (j.lU) 
would  not  "be  used.   They  must  "be  replaced  "by  more  general  solutions  of  equations   .  . 
(7.10)  wherein  S  euidyU  ajid  have  functional  forms  determined  "by  the  non-homogeneity 
in  the  layer. 

Of  course, all  of  the  theory  of  this  article  applies  vhen   the  source  is 
a  dipole  and, in  fact,  simplifies  considera"bly  as  article  7-1  illustrates. 

Ifh=Oorif— <<lwe  have 
A 

©2C0)  =  0 

Q^Co)     =     CO 

and  thus  o"btain  from  equation  (7.S)  the  simplified  expressions 

'  z 


2 


^  (7.15) 


V/(x,y,z)  =  0 

Thus  equation  (7.I5)  with  equations  (5*3)  gives  the  field  of  an  ar"bitrary  source 
situated  011  a  perfectly  conducting  (flat)  earth  above  which  there  is  a  nonhomogeneous 
stratified  medium. 

7.3  The  Special  Case  of  a  Multilayer.  \ 

In  the  event  that  either  or  "both  stratified  nonhomogeneous  media  consists 
of  a  num"ber  of  homogeneous  layers  separated  "by  parallel  planes  the  above  theory 
applies  since  the  discontinuities  in  £  ,  La.,   and  O"  are  of  the  finite  jump  type,   V/e 
shall  apply  the  "basic  theory  of  article  7  "to  this  case. 
:.,v^y..-         To  determine  the  field  vmich  arises  from  an  arbitrary  source  we  must,  in 
"aCi'pr dance  with  equations  (7.S)  and  .  (7.9 ),  determine  first  the  ©(z)  and  G(z)  above 


and  belov;  the  source.  Let  us  suppose  th^t  the  source  lies  below  the  layers  so  that 
the  effect  of  the  la; 
theory  of  article  f , 


the  effect  of  the  layers  will  be  on  the  form  of  0-,  (z)  and  0,  (z)  in  the  notation  and 


37. 
In  accordance  with  our  general  theory  v/e  Eust  first  solve  the  Eiccati 
eqiiation  (7, 10),  namely, 

A  — 

in  the  region  z   >  0.      Since    £  ,   a<.  ,   and    C  nay  have  finite  discontinuities  as  ue 
pass   from  one   layer   to   another  it   is   necessary  to   olDtain  a  separate    solution   for    O 
in  each  layer, though  as   remarked  in  article  6,G(z)  \/ill  be  continuous  throughout 
the  nraltilayer. 

Let  &^,  /^  ^  "ce  the  material  constan-cs  of  the    V-th  layer  and  let  X^ 
be   its   thickness.      Since    ^^^and  Vvi  s-"<3   constant   v/ithin  the    V-th  layer  v;e  can  solve 
for  ©directly,   as   follows.      "Je  first   restate  equation   (7. 10)   thus: 


0' 


r 


2ni 


^2  X     '^v 


1-  -I©' 


£, 


1-  ^  © 

1  *    -^  0  ^ 


Y 


gni 


=  Ce 


18 
v/here  C    is   a  constant   of    integration      ,      If  v;e   compute   the   definite   integral   from 

(3  .to  0v.^.-i  .    ^^s  values   of  O  at  the  bottou  and  top  of  the    V-th  layer,   ve  obtain 


V^l^  (7.17) 


^^ 

'  - 1  '^^ 

C  is  not  arbitrary,   "ve  do  not  require  its  value  in  the  steps  iircnediately  follow- 
ing but  it  can  be  determined  by  the  condition  that  ©(z)  in  the  V-th  la;/er  must  be 
continuous  v/ith0(z)  in  the  (  \/+l)-st  layer,   Finally,  as  the  next  few  steps  indi- 
cate, 0(z)  in  the  uppermost  layer  is  constant  and  determined. 


By  solving  for  c?^  and  using  trigonoinet,ric   identities,   v;e   ottain 


©  ^,    -.  i  -1:^  tan  ^   <  i 
V 


1  -  i 


e^  = ): (7.18) 


which  allov;3  us  to   determine   succecsively  the  quantities  O  i">oa  the   last  one 

'"uv.iich  is   the   solution  of  the  "basic  equation  (7.10)  for  ©in  the  semi-infinite  region 

zyil-.'^y^'^       "^  A    ,   ^^herein  Y"  arid  ^5  have  the  values    V'    and  £.    . 
12  HI  mm 

A  similar  formula  is  obtained  for     9^,   namely 


iti^^- 


0 

P  ,    ,    -  i      ^^'-     tan     f   V,  i^ 

/^m 

(7.19) 


We  can  now  use  these  fonaulas    (7. IS)  and  (7»19)  to   obtain  some  of  the 
quantities   required  by  the  general  formulas   (7.S)  and  (7.9)    0^  t>y  any  specialization 
of  these  formulas.      These  fjeneral  formulas  call   for   the  values   O^(o)  and  ©,(0).      In 

oup  ease,  where  the  multilayei"  lies  immediately  ^bove  the   source,  <3,(o)  =   0     and 

I  0 

0- (0)  =  G  .   The  quantities  0  and  C  can  be  obtained  by  computation  with  the  re- 
cursion formulas  (7.12)  and  (7..19)- 

The  general  formulas  {J.S)   and  (7.9)  require  also  the  functional  ex- 
pressions 0^Cz)  and  9^(z),  O  (z)  ?nd  G  (2).   In  our  case,\;here  the  mxiltilayer  lies 

immediately  above  the  source,  0p(s)  and  QAz)   are  unopecialized  and  depend  upon  the 

medixim  belov;,  ©,  (z)  and  0,  (z),  however,  depend  upon  the  multilayer  and  can  be 

obtained  from  formula  (7.I6)  and  the  corresponding  formula  wherein  /^y  replaces   €.  . 
(see  article  ],h).     \le   note  that  the  expressions  for  0,  (z)  and  ©-.(z)  \;ill  vary  from 

loyer  %o   layer  because  ^x/tM^ft   a^d.  Yi.   ehcinge.  As  a  consequence  the  expressions 
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>v>for  V  and  W  given  "by  formulas  (7.8)  and  (7.5)  will  change  from  layer  to  layer  , and 
"'"•;ii,„t'inately  so  will  the  expressions  for  the  components  of  E  and  H  as  determined  "by 
equations  (5.3)» 

7,U  The  Dipole  in  a  M\iltilayer  of  Three  Strata 

The  formulas  of  article  7*3  ?illou  us  in  principle  to  express  the  solu- 
tion of  the  case  of  the  juultilayer  by  a  series  of  integrals.  However,  the  intet^rals 
"become  complicated  rapidly  if  the  number  of  layers  is  large,   j.'e  treat  in  the  follow- 
ing a  simple  example,  which  represents  roughly  to  a  certain  extent  the  situation  in 

the  atmosphere  of  the  earth.  Consider  a  medium  of  three  strata 
2 

or  a)   0  ^  z  <  i:  ^  ;  £^  fixed  and  ^^  =  1 


-g 


^^J_  b)  i?  i  z=    ^o+K'h   ^i^ed  and/<i  =  1 


^  OJ  <o 


The  ground  surface  is  assumed  to  be  perfectly  conducting  and  we  shall  therefore  use 
results  of  article  7*2,   ^^le  source  shall  be  a  dipole  oscillating  in  the  direction 
M'  the  z-axis  and  located  at  (o,o,o)  so  that,  in  view  of  formulas  (5*7)  ''^^^   (5*2), 


V  So  ^     ],   rr  ik(c;:i:  +^0y  -^   z  /^^  -^  ~ {^^^   ) 


Vi   =  0 

if  the  medium  in  which  the  dipole  lies  is  homogeneous  \/ith  <£=  £  « 

3y  means  of  formula  (7.15)  '^^   obtain  the  solution  for  the  above  inhomo~ 
geneous  mediuo* 


I """^2 2 

Since  V  -/£■  LL     -  c!L     ~  (h        »  '^'^^   expression  for  V  becomes 

r 

^         ff    .        ik(<xx  +By  -^I £  O  (z)  dz) 

V(x.,y.z)  =     frX~  //      V^  '^^^^  ^'^'^^'^ 

The  function  0(z)  v/ill  be   foiind  by  the  aid  of  the   formula-  (7.16). 


IrO. 


Before   doi;ig   this  let  us   coriSider  the  special  case  of  a  corapletely 
honogeneous  uediua  c=     £.     a"bove  the  conductinc^  ground  surface.      Me  have 


and  hence 


V(x,y,z)  =    -~ 


e 


y£o-=<'-^' 


=    2 


e 

ikr 


i.e.    the  original   source  \;ith  four  tines  the  iiitensity     (  lv\    )   . 

We  return  to    the  case  of  the   three  layers.      Since  \ri. thin  any  one  layer 


Porziula  (7.20)  nay  then  "be  rewritten  and  transfomed  as  follows:      'ie  have  directly 

1      f  f   1       ik(c.x*/3y)     ik^'c5  6)(Odz 
V(x.y.z)  =     ^  I   I'q^^  '  ^  dc?4    dp 

We  noi;  make  the   following  changes  of  variables. 
Let 

c><    =    ^   COS   UJ  ^ 

/3    =    /=*  sin  Cf> 

X     =     r  cos    \ij 

y  =  r  sin     "u/ 

The  change  from  CA  ,  /^  to   P  ,    ■J    involves  a  change  in  the  vaxiahles   of  integration  as 
was  performed  in  article   5.     The  change   froia  x,y  to   r,  'XfJ does  not    involve  the  vari- 
ables   of   integration,  ■    •' 

V/e  obtain  .    ,^ 

2n  03 

dz 

/^d  P  d  (p 


Ul. 


"■".'■e  nov;  use   the   fact   that  0(z)  depends   only  upon  Pto  write  our  integral  thus; 

1-:/    £0(z)   dz 


00 

r  - 


,Tr 


V(x,y,2)  = 


TT  e 


} 


WF 


^ikrPcos  n       ^  P^LP 


-TT 


v/herein  n  =  W-   W'and  \/e  iiave  ch?jigcd  the  bounds  of  integration  of   n  to  -n  to  tt, 
'./hich  wc  nay  do  "because   the   integrand  is  periodic. 

1° 
V/e  nov;  use  a  standard  integral  expression  for  J   (x) ,    nauely, 

T   ,     \  1         /  '^'■^  COS  1:      ,    ,  , 

-TT 

to  rewrite  our  integral  as 


V(x,y,2)  = 


CO 


J  (k  r   r)  e 


ii^/  s.e  (2) 


dz 


p±'^ 


(7.21) 


To  evaluate  V  we  i.-nist  first  detemune  O  (  P  )  auad  the  value  of  the 

/  ^©Cz)  dz.  We  proceed  to  do  this»  For  the  determination  of  (^     we  have, by  the 
recursion  fornalas  (7. IS),. 


o,  -  -7^  =  A  -  r' 


Qi  = 


TT  A  '""  F  ^1  ^'1 


1  -  i     ->- 


/. 


^n  41    ^.    X^ 


A 


1    '^  1 


ai- 


i  A5.     tan     ^    -^      ^7 
^  1^     ^"^     A        -^o    ^  0 


©0  = 


i_  i  -i2.  a  tan  in  v;e 


^     -'I    '^    X 

2o 


o     0 


bmce 


v^i  =  /^i-r  -'  ^  ^1- 


1>   s 


1  +  s 


•^^'Stratton,   loc,   cit. ,   p.    ^Sj . 


(7.22) 


U2.        . 
and  since  only  the  even  functions   of  either  radical,    na^nely, 

enter  into   forijulas    (7.22),    the  value   of    0   (d)   for  the  positive  or  negative  value 

of  the  radical  rust  te  the  sane.      Hence  0      is  a  one-valued  function  of   s,    i.e. 

o  ' 

tlie  Rienann  surface  of  Q    (s)  has  only  one  sheet.   The  only  possitle  singularities 

■  0  ■ 

of  the   function    — =- — |i — r     thus  are  poles.      Their  location  is  given  by  the  roots  of 
the  equation  Q    (s)  =  0,    i.e.   "by 

i        »/     ^  ?T\   ^      ft 

.00.  2tt-/'<5_  ^1 

i  tan    -^  ^0   ^0  -     T- : 

1  _ ,  Yi '  '^'*  ~5r  ^1  *^i 

\\q   next  determine  the  function 

ik^^'cGCz)  dz 

v(z)  =  ~-    e  (7.23) 

^0 

in  formula   (7.21).      V/'e  can  do  this    oy  first   solving  the  Riccati   equation  (7^10)   for 

O  (z)  and  then  carrying  out  the  integration,.  Hov;ever,_  it,  is  simpler  to  proceed  as 
follo\7s.  Ve  know  oy  fomijla  (t).5)  that  v(z).  must  be  a  solution  of  the  differenhial 
equation 


t    \ 


l{-^)      ■^k'^   (£    -    /?'')  V  =  0 

Since    £(z)    is  constant    in  tne  individual   layers,   uc  may  treat   this  differential 
equation  as   an  ordinairy  second-border  equation, and  we  conclude  that  v(z)  must  have 

the  forn 

ik    ^2  -ik    S  z  , 

.-A  0,^  O  <  A 

1     -     A^  e  +  ^^   e  ,.  0  =  z   <  z^  =   j>^ 


ik  ^,(z-z    )    ,  -ik    t -X-L-i^)       .,  ^  o  . 

i  ,    .        .  (7.2U) 


ik  ^  o(z-z^) 


i^  =  A2  e  -         2'  .    z^  i     z 


^3. 


D 


<i  •;, 


I, 


JL 


The   special  ior..i   of  v(z)    in   the  nedi-ur.!   z    >  z      is   a  con.- 

seQuence   of  the   ■boundary   condition  that  ©(z)   -7>     5 /£ 

1         V  (  2  ) 
as   z  -»   CO.      Since  G(z)   =  -^      'l^\z)        ''"'^^   condition 

v'  ■/ 

IS     1_     _   ik  -|_  V  ->  0   if   z   ->   00.      iTou  equation   (7.23) 

it  follo-./s   that 


v'(o)  =  ik£ 


0  - 


and  hence,    fron  the  first   of  erpressions    (7.2'4), 

A.  -  B„  =     ^  (7.25) 


1 
V;e  now  inpose  the   conditions   (see  state:.ient   (5«l))   that  v(z)  and     —  v'(z)  i.iu£t  te 

continuous   every\-/here.      These  conditions  require  that  the   first   t\/o  expressions   for 
v(z)  and  v'(z)  in  (7.2U)   r^ast  a^iree  at   z^   =    x    ,and  the  latter  tvm  expressions  and 

their  derivatives  ntust  agree  at     z_  =     x      +    .(,  .     This  gives  us    a  systeu  of  four 

linear  equations   for  the   quantities  A   ,  B   ,  A,  ,   3,  ,    and  A-,   naniely, 


A     e  °       °     4-     3     e  °     °     =  A,    ^  3, 

0  0  11 


ik  i     I  -ik  /     i 

^1       ^0    ^ 


~  -  c         -^ 

0        0        ^_  O'^O  <--o''l 


(A.    - 


^1^ 


A,  a 


ik  -jr^L  i^ 


-ik  /,i, 


(7.26) 


Al   ^ 


ik  ^,^1        -ik2r,i^.,    £^  ^, 


-"l    '^ 


<:         •/         2 


or 


-ik    -C       ^     \  '   c  /, 

OS        -^  O'-Oj.  to  «lv         . 

2A^  -   e  J  (1  .  )     ^^^ 

"- 1  "  0 


•^1       ^0  J 


23     =   e 


ik  V  ^.    l^    I 


0  0 


(1  - 


Sr.         V. 


-)     A,    +   (1 


^o      A 


)  3, 


1 fl       V/       -1  -  £,        r/^l 


1         »0 


2A,    =   (1 


^1     ?^2,    ,        -i^^-    ^^^1 


6i      X'p  ik   /      i 


Those   eouations  allovr  us    to   ejrpress   the    quantities  A    ,   3    ,   A-,  ,    and  5,    in  terr.-.s    of 

0*       0         1'  1 

quantity  A^    : 

A^  =  57    A     e  ;  (1  +  -J-    —r-)    (1  +     -7- -:7-)   e 

-1     Co  ^2       3^1  J 

B.  =     |-    A     e  J  (1  -     — --)    (1  +  -7 7-)   e 

^2  fl     a^  ^2    ^1 


0 


^1     iio  '-2    ^1  J 


(7.27) 


A,=     1    A     e"''^^^^l.'    ^    ^) 
^  2     ^2   *^  ^  £2       ^1 


3_  ik  ^Ti^i  li     / 

B     =      2     ^2   °  (1   -     -^    -^) 

^  <^       <-  ^  2       ''  1 

By  subtracting  the   second  equation  froi.i  the   first,    replacing  A     -  3      "by 

^0  /  \ 

— T—     in  accordance  with  for:.iuln    (7.25)  and  solving  for  A„  ue  obtain 

^0  ^  2 

,      -        1 

where    A  is    defined  "by  the    equation 

\  i^Tc     =    (1  .     l£    ^h    .X  .     Al  12)   e"''^    '^°^°  ^  ^^^^^ 
^0  ^  £1       ^0      '  ^-2    ^1 

^'       ^°  ^'     "^^     '  (7   28) 

-   CI  -  _  -— )(l   +  _ _)   e 

^1     2^0  ^2    ^1 

<^1     tfo  ^2    ^1  • 


r 
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By  using  the   fai.iiliar  relatione   "betueen  trigononetric  and  exponential 
functions   one  verifies   readily  that  5c  can  "be  v;ritten  as  follows,    letting 


^P=  Vi-  r 


£2  =  1  ; 


S2  = 


/o^l 


COS  k    ^q£q  cos  k   2'^^^  ■*•     -zr-  —    sin  k   '^ ^t ^  sin  ^  ^^^^^  j         (7 

0  1  ^^  0 


.2S) 


-i 


—2.     sin  k  ^^i^  cos  k  ^,  i     + 

P  00  X       i 


/l 


cos  k 


i»     ^      sin  y.   0  ^  £ 


0  '^o 


1^1 


which   shov/s    that    ic  (s)    is   an  integral    lunction   of   s,    i.e,    regular   in  the  \;hole 
s-plane,      With  the  aid  of  this   result   \;e  uay  finally  return  to  fornulas    (7,2U)  and 
express   the   function  v(z)   itself   in  the   follov/ing  explicit  fom: 

for  0   <  2   <  E.     =   1    : 
i  0 


v(z)   = 


5^  ( s )  <    cos  k    ^^(z-z^)  cos   k  "^f  ^if^+ ^ -^  sin  k    j/^(z-2^)   sink   ^-^^-^ 

L  1    '^0 

I 


^1  /)  B 

—  sin  k  If  ^i^  cos  k    /^(z-z^)  -     -2.  cos   k  sT^i^^sin  k  2^^(2-2^) 


(7.30) 


for  2^   <   z   <  z 


v(z)    -       7=77 r  COS 

UK  3) 


for  z  >   z. 


k    /^(z-z„)   -^   i   s     -r-  sin  k    ■o,(z-2_) 
1  (i  Si  -^  ^ 


^  iks(z-z    ) 


All  these  functions  are  one-valued  functions  of  s  =    yl  -    P   .     The   only 
singularities  are  poles,    given  by  the  roots   of  the   transcendental    equation 


hG. 


The  function  V(x,y,z)  itself   is   given  Tsy   the  integral    (7.21).      V/e   introduce   the 

variable   s  =  >Jl  -  P       instead  of  p  and  obtain  the   solution  of  our  problea  by  the 
integral 


/■ 


VCx.y.z)    =     -pr-l    v(z,s)      J    (IZT  Jl   -   S   )l 
0  iJ 


ds 


Mhero  v(z,s),    defined  in  fornula   (7.23)    is  given  by  formula   (7.30)'      ^^e  integrand 
is  a  one-valued  and  neroinorphic   function  of  s  with  poles  as   singularities   only. 
The  path  of   integration  in  the  complex  s-plane   is   shov/n  in  the  figure.     Any  deforma- 
A^aT /~t  L.  tion  of  the  path  mst   consider  the  possible  poles   of   the 

Soi~/£>  •>  integrand.      Of  special   interest  are  the  values   of  V(x,y,z) 

B.%    the  boundary  surfaces.      Here  we  have  to   ir^sert   in  our 
integral  the   fuixtions 


z  =  0;    v(o,s)   =        -  /-•■Y  i  cos  k    6  „  A   ^°s  k    {     i     -^  ~  -—  sin  k    /^  y.      sin  k  ^   i 
jiUi  o'Co  1      1       £^    ^^  00  11 


.^1      . 


-is(—  sin  k   /^X-,_   cos  k    ^^i:^  +  "^Aios  k    /^^^  sin  k  /  ^/^)1 
=   z^;    v(z^,s)  =     -^YTiT  ^^°"  ''^    A^l  ~  ^^     ~  ^^"  ^    "^''I'^l^ 


^^  ^2-  -("^2'^^=    -TTTiT 


'■•'^p.  .particular  ,\;e   obtain  the  function  V(x,y,z)  at   the  upper  boundary  by 

'■>■■'■    '"1 


V(x,y,z    )  =     i—    ;  ,    ,     J   (kr  >/l  -  s^)    s  d  s 

The  electronagnetic  field,   finally,    is   found  by  applying   the  formxilas    (5.3) 
v;ith  W  =   9   since  our   source   is   a  dipole,- 


H,    =  -ikV 
^2  =   ^^^x 


tS,   =-(V       +  Y      ) 
3  XX         yy 


Hj^O 


^7- 


With  our  results  \/e  r.av  e:-'-press  the  solution  for  the  sicpler  case  uhere 
;he  aediun  consists  only  of  tuo  strata  separated  "by  the  plane  z  =  z  .  Ue  obtain 
Z-  the  solution  ty   assuLiing  K -^   =  0,  i.e.,  "by  introducing 

in  the  integral 


A, 


e 


V(x,y,z)   = 


the   functions 


So^: 


•(z,s)  J   (kr  >/l  -  s^    ) 


S     OS 


Ho 


<   z    <   z^;    v(z,s)   =      -3:^^     (cos   k    Yq(^-z^)   ^^   ^^   TT     ^"-^  ^    ^c^"~^l^ 


^1  <  ^'  ^^^»'^  =  -jhji 

whereS^Cs)   =   s   cos  k  'Y 
i.e. 


^  ■. 


V(x,y,z) 


cos   k  V   (z-z-,  )   +  i 


ik  s(z-z,  ) 

sin  k    V^(z-Zt  ) 


(7-31) 


-  1  — r—  sxn  k    *      ..  ^ 
0  '^  o  C  0  ''  0 

*-o 


^0 


L  s    cos  k  Y      X     -  i 

o      0 


^  sin  k    Y    J^ 
£0  °     ° 


JgCkr  v/l-s    )s   ds 
(7.32) 


, Tj-     ik   s(z-z    ) 

J  (kr /i  -  s"^  d-         "■  ; 


s   cos 


!•:   Y.  i? 


O  0 


s   ds 


0     0 


sin  k   V     X 
0  '  0 


\/here  the  first  expression  is  valid  for  0   <  z   <   z,  ,    the   second  for   z   >  z,  . 

Articles  ] ,h  g?    j' .  7  gg'^.TT^.-'f^    tliat  the   dipole   \.'as  on  a  perfectly  conducting 
ground  and  therefore  utilized  equation   (7.I5).      Hox/ever,    in  vic\;  of  the  discussion 
preceding  this  formula  it  \/oald  be   equally  possible   to   consider  the    case  of  the 
dipole  above  any  sort   of  ground. 
7.5  £xaaple    of  a  Stratified  i^iediuia  v/ith  Continaously  Varying  Index   of  Refraction. 

Ihe  preceding   article   showed  ho;;   the  general   theory   of  an  arbitrary   source 
located  bet\;een  tv/o   semi-infinite   stratified  media  could  be   specialized  to   the   case   of 
a  dipole  located  on  a  perfectly  conducting  ground  with   the  medium  above  ground  con- 
sisting of  parallel  layars   each  v;ith  constant    £  ,  /^  ,    and    C5  .     The  present  article 
v/ill  indicate  hov;  the  general  theory  can  be  specialized  to  the  case   of  a  continuously 
varying   dielectric   constant.      For   the  purposes   of   illustration  v/e  consider   the   strati- 
fied medium   in  which yU  =  1  and  £.  varies  according  to   the  expression 


£  =  1  + 


k^z^ 


(7.33) 


us. 


v/here  p  is   a  positive  or  negative  constant.     V/e  assume   t}iat   our  medi-um  is  "bounded 
at  z  =    O     >  0  by  a  homogeneous  medium      £,-,      A^-,.so   ti:iat  the   singularity  of   forcula 

(7«33)  at   z  =  0  is  excluded.      Since    Q    and  p  are  arbitrary  constanta  ^;e  can  approxi- 
mate any  irJaomogeneous  medium  in  uhich  c  increases   or  de- 
creases monotonically  "by  the   formula  {7»35)  with  the  aid  of 
the  tuo  paraxaeters  p  and  o  .      7or  example,    if  the  values 
£,    and    £,p  at  the  hottom  z  =6   and  at    the  height   z  =  H  -t-o 

are  kno\m,ve   can  represent   the  medium  hy  formula  (7»33)» 
usin^r  the  parameters 


o 


->x 


s  = 


-  1 


-  1 


(7.3U) 


p  = 


(  £^  -  1)  H^  k^ 


-  1 


\Ie   assume  that  a  source  is  located  at  z  =  h  >  0   \/hose  radiation  pattern 
is  kno\m  if  the  surrounding  medium  is  homogeneo\is .   In  particular  v;e  shalD.  assume 
that  the  function  V  of  the  two  characteristic  functions  Y   and  U  is  identically 
zero,  so  thatthe  electromagnetic  field  has  the  special  form 


\ 


=  ik  W 


=  -ik  v; 


h 


=  w 


sz 


X 


^3   ° 


2    yz 


3     yjc   yy 


(7.35) 


It   is   clear   that   this   is    a  very  special   class   of  transmitters,     IVc   choose  this   class 
"because  the   differential   eo^uation  for  '.i  admits   simple   solutions  under  the  assumption 

(7.33). 

According  to   formula  (5.7)  "©  thus  consider   sources  which  have,    in  a 
homogeneous  medium,   the   radiation  pattern  Y  =  0  and 


for  z  >  h:     '■/,   = 


for   z   <  h:      V/     = 


do.    d/i    . 


jl    w^(<^.^,z)e^^^^^^*/^^) 
/Tw2(cX,^,,z)   e''^(^^^--,^y^      ac    dr3 


(7.36) 


vmere 


w^  =  G^(  a,  B  )e 


^9. 
ik(z-h)    s/n^-  c^^  ~    p^ 


& 


w.  =  &^(  a  ,    3  )e 


v/here   " 


n     =     £      /x     ,    £,       and    £4      Iseir.g  characteristic   of  the  honogenecus  medivun. 

G,    and  Gp  being  given  functions,   the  modifications   of  forcrala   (5«7)   ^.re  required 
only  because  the  source   is  at  z  =  h,      ?roir.  their  very  definitions    it  follous  that 
\t     and  w     satisfy  at   z  =  h  the   follo\;ing  characteristic  discontinuity  conditions; 


w^(h)  -  v/2(ii)  =  \  -  ^2 


(7.37) 


w|(h)  -  u^(h)  =  ik,/n^  ~  p^   (G^  ^  G^) 
By  these  conditions,   and  by  the   conditions 


-7—    v;,'(z)  -  i  — —    v/n^     -P         \:.'\z)  — >  0     if  z  -^  co 

To       ^  Mo    ^        ^  '  (70S) 


~ —    w'Cz)  +  i     —    Jn^  -    p^     vJz)  -^if   z  ->  - 


00 
0        2  /"^  0 


/..        . 


which  follow  from  the  results   (7.11)  v/here  9(z)   =  ~-rr-    — —7    >        (cf.    formula 

ik        /.^  \/   ( 2  ) 

(6.U)),and  \;hich  must  hold  also  when   the  mediuEi  is  homogeneous,    the  solutions 

v-(z)  an.d  w   (z)   of   the  ordinary  differential  equation    (See  formula   (5«5)) 

(liL_)'.k2(n2.^2)v=0  (7.39) 


are  uniquely  determined,  if  the  medium  is  homogeneous,  Ue  thus  may  characterize 
the  radiation  pattern  of  our  source  by  the  discontinuity  conditions  (7.37)  and  the 
boundary  conditions  (7.3S). 


50.- 

Had  v;e  considered  the  ciacs  of  trancmioterc  for  \;hich  V/  =  C  <'^nd  V  given 
"by  an  expression  such  as  (7*3^) »  ^  class  which  includes  the  vertical  dipole,  v/c 
should  "be  led  to  an  equation  like  (7*39)  except  that  y.i  \;ould  be  replaced  by  t   . 
Since, hov/cver,£  is  variable  uhereaz  jU  iz   constant,  the  solution  of  the  annlo^-^c  to 
equation  (7.39)  becomes  sruch  more  difficult.   This  case  ins  net  been  treated  as 
yet. 

If, now, the  .viediiiEi  is  inhbmogeneous,  as  in  fornula  (7o3)»^"''3  obtain  the 
electronagnetic  field  by  the  formulas  (7.3^)  "here,  however,  v;,  (z)  and  v;^(z)  are 

solutions  of  the  differential  eo^uation  (5«5)|  namely, 

w"  +  k^  (  £  -  p^)   w  =  0  (7.^0) 

in  which  S.  is  no  longer  constant  but 

e  =  1  +  —P 

,22 

k      z 

and  wherein  we  have  taken  JLX=  1. 

v/-(z)    and  w„(z)  must   satisfy   the   conditions 


v;^(z)  -  ik  yi  ~  f)^     v/^(z)  -->  0  if  z  ->  CO  (7-^1) 


v;^(  ^    )  +  ik//,  yn|  -    yO^     w^CS  )  =  0     at    z   =§    . 


2 


where  n-,  =     & -,  /JL^   is   determined  by  the  rr^terial   constants^    ,   IJ.      of  the  boundary 

mediuQ  z   <  o     »     ^^-e   first  condition  is   rx-reli'  a  restatement  of  the  boundary  con- 
dition  (7,11)    in  terms   of  v/-,(z)    instead  of     Q^{z).      The  second  condition  exioresses 
the   fact   that   the   S(z)    in  the   region   o  <  2   <  h  must  be   continuous   vri.th  the  G(z)    in 

the  semi-infinite   region   z   <  o   .     How, in  the  latter.   0(z)  =  0„(z)   = =^  .  by  one 

/X3  ^^ 

of  the   conditions    (7*11 ).      Hence   the   former  ©(z)  =        ..  — ^^ — {■ — r-    must   equal 2. 

ik^w(z)        ^       it 

at  z  =  0  .   la  particular  w  (  O  )  =  0,if  the  medium  z  <  ^  is  a  perfect  conductor 
v/ith  E-j  =  00.   (See  article  7.2). 

The  conditions  (7.37)  remain  unchanged  in  our  medium;   This  expresses 
the  fact  that  the  same  source  is  placed  in  our  inhomogeneous  medium. 

Our  problem  is  to  determine  two  solutions  v;,  and  w  of  equation  (7.^0), 
v/hich  satisfy  the  boundary  conditions  (7.^1)  and  the  discontinuity  conditions 
(7*37  )•   ^or  simplicity  let  us  assiime  that  the  sourde  is  symmetrical,  so  that  C 


51. 

Cr^  =  a^  =  Qi  ex  ,  f:i  )  =  &{  p  )   v/here,  as  usual,  p  ^  =  d^^  +    jl? ,      Let  us  also  aosune 
that  the  aedium  z  <  ^  is  a  perfect  conductor.  Then  our  problen  is  to  find  t;/o 
solutions  w,  (z)  and  mAz)    of  the  equation  (7.U0),  nanely, 


v-'(z)  *  y?  (1  +    -^   -  p^)  V  =  0       (7.U2) 

k  z 


whidi  satisfy  the  conditions 


i/^(z)  -  ik  Jl  ~    p'^     v^{z)   -A  0,    if  z  -^ 


CO 


v^io    )  =  0 


(7.U3) 


and  the  conditions  at  z  =  h 


w^(h)  -  v.'^Ch)  =  C 


„^(h)  -  u^Ch)  =  2  ik  rV^Zp^  G(  p  > 


il.hh) 


where 

r  =  1  +  _P_ 

°        1^2,2 
k  h 

The  general   solutions  of  equation  (7«^2)  are  given  by  the  Sessel 

functions  and  Eankel  functions 


JT-    J^    (k  Jl  ~p^     z)   and    Jir    H,^     (k  71   - p^     z) 
as   one  easily  verifies.      The   index,   v  ,    is   determined  by  the  relation 


p  =  1/U  ~  V^    \ 


I.e. 


V  = :  fTj^rr;  (7.^5) 

We  satisfy  the  conditions  (7. ^3)  by  choosing 


w^(z)  =  0^/7"  ^  h;  (k  z  Jl  -  p^     )  (7.U6) 

/'..2     / 2  5^  X 

v;  (z)  =  Cg/r-(''^  \^   E^^(kyi-£>2   ,)  _h2  (k/l-O^   z ) ) 


where  C^  and  C„  are  arbitrary  constants, 


52. 

These  constants  will  "be  determined  "by  the  conditions  (7.^^).  '''c  intro- 
duce temporarily  the  notation 


5l^(z)  =  H^^  (k  zv/l  -^2  ) 


v2 


i^Sjl-    P^   )     _., 


Then  to  meet   the   first  condition  of   (7.^^)   "e  vrito  ",  (z)  and  ^'Az)   in  the  form 


v;^(z)  =  cy2h    Sl^(z)A.2(h) 

v^Cz)  =  cyih   S^3_(h)Sl2(z) 
where  C   is   still  arbitrary.     The  second  condition  of   (7.^^) 

w^(h)  -  w^(h)  =  C  h^^^(h)rt2(h)  -n2Cii)A^(h))... 


=  2  ik  Je^  -  f?-       G(  f> ) 
then  yields   C.     With  the   aid  of  the   identity^*^ 


^(KtCx)|j     h5(x)-H^Cx)^    Ki^:0)     =-     % 


applied  to  the  first  of  equations  (7.^)  we  obtain 
-Ui 


I.e. 


and  hence 


VC=  2  ik  j£^-f>^     Ci^) 


W^(z)  =  -   n^y£^_^2   ^^^^  /I^9.^{z)Sl^ih) 


"2(0=  ~^y/a^-p^    G{f=)/^Sl^{h)Sl^{z) 


(7,^7) 


(7.^S) 


(7.H9) 


20. 
_^^^^  See  "./atson,  G.  1>I, :  Theory  of  Bessel  Functions,  p.  76,  formula  (U), 


53. 


The  case  Gi  P  )  =     — —  —     is   of  special    interest  since,    for  a  honogencous 

medium   £,=     S_  ^^^ /^q  ~   ^i    ^^^  eo^uations    (7.3°)  yield  the  function   (see  fomula  5.5) 


ikr 

Such  a  source  thus  gives  in  our  irjiomogeneous  mediua  the  functions 

k 


w^(z)  =  -|    /In  ^^(z)P.2C^) 


(7.50) 


w^, 


,(z)  =  -  I    ,/ih  ri2^z)fli(h) 
or,  explicitly,   in  viev/  of  the  meaning  ofic,(z)   andyt-jCz) 


1  H 


/ih  fn^^i^y^  _^^2  ^  _p(^)  H^^  (j^y^  _p2  )y  n^d^^Vi  -^2_j 


(7.51) 


WgCz)  =  ^    y^    f  H^    (kz  Jl  -  p^    )   -   n  O  )     hJ  (kz  Vl   -  /i2      )  j    jji^  (kh  yirrp   ) 

V 


v/he  re 


The  function  v;(x,y,z)    itaelf  finally  follo\;s  from 


\-U  v^.f»«"='°^^^'^^>i^i^  (7.5a) 


for  z  >  h  and  from 


■^'2  =  //    "2(^'/'^    e^^^^^-'l^y)      doCdj3 


for  z   <  h. 


5U. 

The  argument   for  this   result   is   precisely  the  arguiacnt   for   formula    (6.I3).      However, 
in  this   article   the   theory   deals    directly  \;itn  u(z)    instead  of  replacing  an  cx- 
I^rcGSion  in  w  "by  0(z).      From  formula    (7.52)   the   electromagnetic  field  can  "be   deter- 
mined "by  iiieans   of  the   relations    (7-35)« 

Our  result    occoincs   quite   simple   if©  — >  0.      This   limit   will   "be   ootained 
if  ue  assume   that    V   >  0,    i.e.,    that  p   <  l/U,   for  only  in  case  'V  >  0  is  a  definite 
limit   assured  for  /    iP),    v/hich  alone   contains   the    quantity     O  . 

In  order  to  find  this  limit   the  relations   which  follov;  directly  from 

1  2 

the  definition  of  H  ,  and  E^    ,    ncimely. 


2  J        (x)  =  e  h     (x)  +  e  h  .  (x) 

—  \/  V  V 


2  J^  (x)  =  hJ,  (x)  +  h;(x) 
are  used.   Since  J  .  approaches  0  as  z  '^  it  follows  that 


(7.53) 


•^v  ^^)     _        1 .  r 


J     (x)  iVTT_^    -iVTT/'^ 

~  \/  e      +  e       / 


The  limit   of  the  left   side  as  x  approaches  0  is    zero    if     V  >  0.      Hence    f^—^     -1. 

Hence  as  O  approaches   0   the   formulas    (7. 51)   thus  "become,    in  view   of  the 


1  2 

definition  of  H       and  H 

V  V       * 


(7.5^) 


"l(^)  =     ^    ^/^^     Jy  (k  h/l  -io2    Jh^^  ik  z  Jl  -f      \ 

"2(2)  =   ~  /^'n   J^U  z  v'l  -^2  Va  ^j,  h^i  -f^   ) 

If  v/e  substitute  these  results  in  formula  (7.52)  and  take  advantage  of 
the  steps  already  demonstrated  in  going  from  formula  (7.2C)  to  formula  (7.21)  -  the 

f 

^  iky  £  ©  (z)  dz 

quantity  w^(z)  here  plays  the  part  of  —^  e   °  there  -  we  obtain 

Oo 


55. 

fe«-.  ■.  b,?  z   >  h 

[  A^^  

\  W^(x,y,z)  =  knjzh  J     (k  h  71  -  /^^  )  H^^  (k  z  ^1  -  p^   )  J^(kr  /^)  f  d  /3 

i 

and  for   z   <  h  (7.55) 

00  

\!^{x,y,z)   =  k-rr  /zhj         J^  (k  z  v/l  -  j^^     )  H^  (k  h  Jl  ~  p^  )  J ^{kr  f^) pd  f> 


I  Z  2 

where     r  =  >/x     •♦■  y       . 

This  result  can  perhaps    oe  further  Gimplified  "by  the  use   of  properties 
of  Bessel  functions.     However,    this  \/ill  not  "be   done  here,  for  the  primary  purpose 
of  this  article  is   to  indicate  hoi;  the  general  theory    specializes  for  the  case   of  a 
medium  with  a  continuously  varying  index  of  refraction   (or  dielectric  constant ).      V/e 
are  not  at   the  moment  interested  so  much  in  ^he  particular  law  of  variation  nor   of 
the  particular  class  of  transmitters    involved  as   in  exploring  the  type   of  integrals 
to  which  the  general   theory  leads   in  particular  cases.      Application  of  the  above    ■ 
theory  to  .investigate  lav;s   of  variation  of  the  C  or  of  .the  index  of  refraction  which 
are  actually  found  in  the   atmosphere  hinges  upon  the   direction  of  future  research  as 
■  discussed  in  the  introduction  to  this  paper. 
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